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Abstract 



The overlap Dirac operator, which satisfies the Ginsparg- Wilson relation, realizes 
exact chiral symmetry on the lattice. It also avoids fermion doubling but its locality 
and smoothness are subtle. In fact, the index theorem on the lattice implies that there 
are certain points where the overlap Dirac operator has discontinuities. Aside from the 
theoretical subtleties, this non-smoothness also raises practical problems in numerical 
simulations, especially in Nf ^ full QCD case. One must carefully calculate the 
lowest eigenvalue of the Dirac operator at each hybrid Monte Carlo step, in order to 
catch sudden jumps of the fermion determinant on the topology boundaries (reflec- 
tion/refraction). The approximation of the sign function, which is a crucial point in 
implementing the overlap Dirac operator, gets worse near the discontinuities. 

A solution may be to concentrate on a fixed topological sector in the full theory. 
It is known that an "admissibility" condition, which suppress small plaquette values, 
preserves topology of the gauge fields and improves the locality of the overlap Dirac 
operator at the same time. In this thesis, we test a gauge action which automatically 
generates "admissible" configurations, as well as (large) negative mass Wilson fermion 
action which would also keep the topology. The quark potential and the topology 
stability are investigated with different lattice sizes and different couplings. Then 
we discuss the effects of these new approaches on the numerical cost of the overlap 
fermions. 

The results of quenched QCD in the e-regime are also presented as an example of the 
lattice studies with fixed topology. Remarkable quark mass and topology dependences 
of meson correlators allow us to determine the fundamental parameters of the effective 
theory, in which the exact chiral symmetry with the Ginsparg Wilson relation plays a 
crucial role. 
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1 Introduction 

Lattice QCD has played an important role in elementary particle physics, in particular, in 
studying the low energy dynamics of hadrons. One can nonperturbatively calculate hadron 
masses or decay constants with Monte Carlo simulations. As compensation, however, the 
lattice discretization of space-time spoils a lot of symmetries of the gauge theory. Violation 
of the translational symmetry would be an easiest example. 

It is well known that chiral symmetry is not compatible with the absence of fermion 
doubling, due to the periodic properties of the lattice Dirac operator in momentum space 
[1,2]. A popular prescription to avoid the appearance of unphysical modes, or doublers, is 
adding a so-called Wilson term [3] to the naive subtraction operator; 



where r is Wilson parameter (we set r = 1.). ^ This term gives a large mass to the doublers 
which are decoupled from the theory. The Wilson term, however, violates chiral symmetry. 
A well known difficulty with Wilson fermion is renormalization. The complicated operator 
mixing, additive quark mass corrections, have to be calculated. Actually these are obstacles 
to obtain reliable numerical data, especially in the chiral limit. 

With the overlap Dirac operator [4,5], as well as the other Dirac operators [6-10], which 
satisfies the Ginsparg- Wilson relation [11], one can construct lattice gauge theories which 
have the exact chiral symmetry [12]. Although the good chiral behaviors in applying the 
overlap operator to QCD are reported in both theoretical and numerical studies, its locality 
properties and smoothness with respect to the gauge fields are not so obvious. Since it has a 



term proportional to = 1/ ^{'^^{aD\Y — 1 — s))^, where s is a fixed parameter in 

the region \s\ < 1, the near zero modes of i^iy can contaminate the locality or smoothness 
properties. In fact, it is not difficult to see that there exist some points where the overlap 
Dirac operator, aD, is not smooth by noting the fact that the index of aD can take integer 
values only. Also practically, near zero modes of Hw causes some problems in the numerical 
^The other notations used here and in the foUowing of this paper are summarized in appendix IXI 
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simulations. Small eigen-modes of \Hw\ lower the convergence of polynomial or rational 
expansion of l/-y/oXff^. For example, to keep a certain accuracy, the order of the Cheby- 
shev polynomial has to be proportional to 1/Amin, where Amin is the minimum eigenvalue of 
In full QCD with the dynamical overlap fermion [13-24], one would have to carefully 
perform the hybrid Monte Carlo [25] updating near H^y ~ points since sudden changes of 
the trajectories, reflection or refraction, should occur due to sudden jumps of the fermion 
determinant. Thus, at least, the smallest eigenvalue of i^vF always needs to be monitored 
in conventional methods of lattice QCD simulations, which is very time consuming and it is 
known that reflection/refraction itself has systematic errors when one employs the pseudo- 
fermion method [23]. Because of these difficulties, no lattice QCD study with the Nf ^ 
overlap Dirac fermions has been done except for the cases with a very small lattice size. 

An interesting solution might be prohibiting the topology change along the simulations. 
It is known that under a smoothness condition on the plaquette variables [26-29], 

||1 — P^,^(x)|| < e for all (x, /i, z/), (1.2) 

which is called the "admissibility" bound, any eigenvalues of Hw are non-zero (we denote 
\Hw\ > 0) and the topological charge can be conserved if e, which is a fixed number, is 
sufficiently small [29,30]. Furthermore, when \Hy\r\ > 0, the locality is also guaranteed. The 
"admissibility" condition, Eq. p.2p . is automatically satisfied if one takes a type of gauge 
action which diverges when ||1 — P^u{x)\ \ — > e [31-38]. 

Topology transitions can also be suppressed by including the factor, detH^, in the 
functional integral [39]. The inclusion of this factor was previously considered in a study 
of domain- wall fermions [40-42], where the aim was to reduce the effects of the finiteness 
of the lattice in the 5th dimension. If any eigenvalue approaches near Hw = along the 
simulation, the determinant, det H^, would give a very small Boltzmann weight, and such a 
trajectory would be rejected. Since det H^r = det{Dw — (1 + s)/a)^ is equivalent to Nf = 2 
Wilson fermion determinant with a negative cutoff-scale mass, it would not have any effects 
on the low energy physics. Moreover, the numerical cost of this determinant is expected to 
be much smaller than that of the dynamical overlap fermions. 

What can we do with the configurations in a fixed topological sector ? A straightforward 
application would be QCD in the so-called e-regime [43-48], where the linear extent of the 
space-time is smaller than the pion Compton wave length L <^ I/itlt^. In this regime (though 
it is an unphysical small- volume situation), it is believed that one would be able to evaluate 
the pion decay constant and the chiral condensates, which are the fundamental parameters 
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of the chiral perturbation theory (ChPT) [49,50]. They should be evaluated in lattice QCD 
studies without taking the large volume limit, since the finite volume effects are already 
involved on the ChPT side [51-66]. 

In this thesis, we study 

1. The practical feasibility of the topology conserving gauge action which keeps the "ad- 
missibility" bound, Eq. ()1.2|) . as well as the Wilson fermion action with a large negative 
mass. A careful analysis on the gluonic quantity and comparison with that with the 
standard plaquette action have to be done. 

2. How much stable the topological charge can be, with these topology conserving actions. 

3. Their effect on the numerical cost of the overlap Dirac operator. 

4. The determination of the low energy constants of quenched chiral perturbation theory 
in the e-regime in a fixed topological sector. This study would be helpful when A'^^ ^ 
simulations with the dynamical overlap fermion are done in the future works. 

We start with the theoretical details on the overlap Dirac operator and topology of the 
lattice gauge fields in Sec. El The technical issues of our numerical study is presented in 
Sec. El To test their practical application, the static quark potential with different couplings 
and e and with/without the negative mass Wilson fermions is investigated (Sec.jH). We study 
the parameter dependence of the topological charge stability in Sec. Then the effects on 
the overlap Dirac operator are discussed in Sec. (HI The numerical result of quenched lattice 
QCD in the e-regime is presented in Sec. [71 Conclusions and discussions are given in Sec. |H1 

The main papers contributed to this thesis are 

• H. Fukaya, S. Hashimoto, T. Hirohashi, K. Ogawa and T. Onogi, "Topology conserv- 
ing gauge action and the overlap -Dirac operator," Phys. Rev. D 73, 014503 (2006) 
|arXiv:hep-lat/0510116| [37], 

• H. Fukaya, S. Hashimoto and K. Ogawa, "Low-lying mode contribution to the quenched 

meson correlators in the epsilon-regime, "Prog. Theor. Phys. 114 (2005) 451 |arXiv:hep-lat/ 0504018' 
[59]. 

Refer also [31,32] which are similar studies in 2-dimensions as a good test ground. 
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2 The overlap Dirac operator and topology 



The overlap Dirac operator [4, 5] is defined by 

^ = - I 1 + 75 = I , a=— aHw = 'j5{aDw -I- s), (2.1) 



which satisfies the Ginsparg- Wilson relation [11] 

75D + D75 = aD^.D, (2.2) 

and 75-hermiticity = 75D75. Here s is a real parameter which satisfies |s| < 1. The 
Dirac operator Eg. 1)2.11) is gauge covariant and has no fermion doubling, as one can see in 
the Fourier transform of aD in the free case ^; 

1 f 1 \ "^^^ 

aD{p) = 1 - (^1 - -a^p^ _ ia-i^p^ U+-a^Y^ P^pA ~ ^al^lP^l + 0{{ap^f), 

Pf, = (2/a) sin(ap^/2), p^ = (1/a) sin(ap^), (2.3) 

where the term o?p^ /2 at ap^ ~ vr in any direction gives large mass, so that the doublers 
are decoupled and the continuum limit can be properly taken [12]. The Ginsparg- Wilson 
relation guarantees that the fermion action 

SF = Y.i^{x)Dij{x), (2.4) 

X 

is exactly invariant under the chiral rotation, even for finite lattice spacings; 

^ ^ e*«75(i-ai?)^^ ^^^ei«75_ (2.5) 

Thus the chiral symmetry at classical level is realized on the lattice. 

In order to see quantum level properties of chiral symmetry, let us consider the eigenstates 
of the operator 075/^ [67,68]. Note that if an eigen-mode of 075/^ has the eigenvalue ±2, it is 
also the eigenstate of aD with the eigenvalue 2, which means that this mode has ± chirality. 
It is easier to see that zero modes of 075/^ can be taken as eigen- modes of 75. Every other 
mode of a'-y^D with eigenvalue A„ in the range < |A„| < 2 is not chiral but has its pair 
with eigenvalue — A„ through the Ginsparg- Wilson relation; 

(075^) I A„) = A„|A„), 
(075^)75(1 -a^/2)|A„) = -A„75(l-«^/2)|A„). (2.6) 



^Here we show s = case for simplicity. 
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These modes have another interesting property, 



(A„|75|A„) = A„/2, (2.7) 

which is again obtained from the Ginsparg- Wilson relation; 

= (A„|{75(75a^) + (75a/^)75 -(75a/^)(75a/^)}|A„) 

= 2A„(A„|75|A.)-A2(A„|A„). (2.8) 

Now it is obvious to see 

= Tr75 = ^ {\n\l5\K) + Yl (An|75|A„)+ ^ {Kll^lK) 

A„=0 0<|A„|<2 |A„|=2 

0<|A„|<2 

= n+ -n_ + ~ N_, (2.9) 

and 

Tr a^,D = ^ A„ + ^ 2 

0<|A„|<2 |A„|=2 

= 2(A^+-A^_) = 2{n_-n+) = -2Q, (2.10) 

where N± denotes the number of A„ = ±2 modes and n± is that of the zero modes with ± 
chirality. It is known that the equation 

- -^tT^ctTsj5aD{x,x) = ■7;^tTce^i,paF^i,{x)Fp„{x) + 0(0^), (2.11) 

can be obtained with the perturbative expansion P^i^^x) = e««^('^Mi'(a^)+o(a))_ Eq. ()2.9p and 
Eq. ()2.1ip show that the Atiyah-Singer index theorem is restored in the continuum limit. 
In this way, the overlap Dirac operator properly establishes the quantum aspects of chiral 
symmetry, the anomaly or the index theorem, on the lattice. 



It is, however, not difficult to see that the overlap Dirac operator is not a smooth function 
of link variables. Consider two gauge configurations, one of which U^{x) = 6"^^°^^^ gives the 
index Q = and another is U^{x) = e^^^t^^^'^ with Q = 1, where A^^{x) denotes Lie algebra of 
SU{?>). One obtains a smooth path which connects U^{x) and Uj^^x), for example, t/*(x) = 
gj(iAi (x)+(i-t)AO (x))^ ^j^jjg Q take integers only. There must be at least one jump from Q = 
to Q = 1, at some t along the path, (0, 1). This jump of Q indicates non-smoothness of the 
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overlap Dirac operator. Hernandez et al. [29] showed that this discontinuity occurs exactly 
when Hw has zero mode. They also proved that if the configuration space is limited such 
that \Hw\ > u holds for a certain positive number m, then the locality is also guaranteed"^. 
To see what happens when an eigenvalue crosses Hw = 0, let us rewrite the overlap Dirac 
operator, 

(^^^-^-'^ ] . (2.12) 

'{aDw - 1 - sYiaDw - I - s) J 

Since every complex eigenvalue A^^' of (aDw — ^ — s) makes a pair with its complex conjugate 
A^^*, only real modes can pass through {aDw — 1 — s) = or Hw = 0. Note that any real 
modes are chiral; [Di^, 75]|A^''^') = 0, and these real modes are the eigenstates of D which 
eigenvalues take or 2/a only, which means that crossing Hw = is always accompanied 
by topology changes. 

In the numerical studies, there is a very crucial advantage of suppressing the small eigen- 
values of \Hw\ ■ In order to implement the overlap Dirac operator, the sign function in 
Eg. (12.11) has to be expressed in the polynomial or rational expansion; 

sgn(x) = = Oo + aix + a2X^ + ■ ■ ■ (polynomial) 



,^2 



or 



{x - hi){x - 1)2) ■ ■ ■ , , . 

= ^0 7 TT ^ (rational), (2.13) 

[x - Ci){x - C2) ■ ■ ■ 

whose errors are controlled to some desired accuracy. The order of polynomial or rational 
function one needs is known to be a monotonous increasing function of the condition number 

^=^, (2.14) 

^min 

where Amin and Amax denote the lowest and the highest eigenvalues respectively. For example, 
the accuracy of the Chebyshev polynomial approximation of sign function, sgn.Q^^^{aHw) , 
with degree A^poiy is empirically known as ~ y4exp(— -BA^'poiy/^), where A and B are con- 
stants [53]. Apparently the approximation of finite order A^poiy would break down when 
Amin = 1/k = 0. In the hybrid Monte Carlo simulation with the dynamical overlap fermion 
(See appendix IHl) [25], the problem would be worse, since one has to consider the disconti- 
nuity of the fermion determinant [16-22]. Every molecular dynamics step in the simulation 
trajectory, one would have to monitor the nearest zero-mode of Hw, and judge if the topol- 
ogy change occurs or not. Then careful recalculation of the link updates has to be done 

''Recently it is discussed that the mobility edge of Hw plays a more important role for the locality 
property. See refs. [69-73]. 
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to determine the simulation trajectory to enter the other topological sector (refraction) or 
go back to the previous sector (reflection). This procedure would require an enormous nu- 
merical cost when one tries the Monte Carlo simulation in a large volume. If one omits 
this step it would make the acceptance rate very low due to overlooking sudden jumps of 
the determinant. Keeping topology and assuring smoothness of the determinant allows us 
to avoid this procedure. Thus, to exclude configurations which gives Hw ~ is useful not 
only theoretically to construct a sound quantum lattice theory with a smooth fermion de- 
terminant (in particular, it is essential for the chiral gauge theory [26-28, 74-86] and also 
applied to the supersymmetric theory [87-90] or non-commutative spaces [91,92] .), but also 
in practically, to reduce the numerical cost of calculating the overlap Dirac operator both 
in quenched and unquenched simulations. Here we would like to make two proposals that 
would avoid the appearance of Hw ~ modes or topology changes. One is the modification 
of the gauge action and the another is the additional fermion action. 

For the former solution, one can construct a gauge action, which generates link variables 
respecting the "admissibility" bound Eq. ()1.2|) [27,33-38]. A simplest example is 

SaJ^%T^jr^^^^F!Mm- whe„ 1 - ReIVP„„W/3 < e. 

[ oo otherwise 

Every "admissible" gauge configuration with very small e keeps \Hw\ > and therefore, 
keeps topological structure of gauge fields. In order to explain how "admissible" gauge fields 
can preserve the topological charge, the best example would be the U(l) gauge theory in 
two- dimensions, for which we can define an exact geometrical definition of the topological 
charge [26,27] 



1 .^1 



2iT V2 



X) 



a'F';Xx) = -t\n{P,,ix)), -n < a' F'^^ {x) < n . (2.16) 

Pfj,u{x) denotes the plaquette in the U(l) theory. In two dimensions, Q^eo gives an integer 
on the lattices with the periodic boundary condition. Since the jump from F'^*(x) = — vr 



to Fj^^{x) = +71 is allowed, the topology can change easily. It is the U(l) version of the 
Liischer's admissibility bound 

1 - ReP^,(a;) < e, (2.17) 

with e < 2, that can prevent these topology changes because the point Fj^^x) = ±tt is 
excluded under this condition. Furthermore, it can be shown that Qgeo is equivalent to 
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the index of the overlap fermion (with s = 0) if e < 1/5 is satisfied. For the non-abehan 
gauge theories in higher dimensions, we do not have the exact geometrical definition of the 
topological charge. But it is quite natural to assume that a similar mechanism concerning 
the compactness of the link variables allows us to preserve the index of the overlap-Dirac 
operator for very small e. Note, however, that the value of e ~ 1/20 is too tight for the 
simulation with the cutoff around 2.0GeV. For practical purposes, we would like to try much 
larger e, as we still expect that the topology is stabilized well in practical sampling of gauge 
configurations, even though the topology change would not be prevented completely. It is 
also notable that the difference between the gauge action Eq. (j2.15|) and the Wilson plaquette 
action is only of 0{a'^) and one can safely take the continuum limit. Also, the positivity 
which would be lost [93] at cut-off scale by the restriction on the configuration space, is 
restored as e/a^ oo. 

For another solution to suppress small eigenvalues of \H\y\, the additional fermion deter- 
minant 

det a^H^ (2.18) 

could be effective [39,40,42]. This determinant would prevent the appearance of near zero- 
modes of Hw, by giving small Boltzmann weights to the configurations which have Hw ~ 
modes. Since it is equivalent to adding 2-flavor Wilson fermion with a large negative mass 
at cut-off scale, —(1 + s)/a, it would not affect the low-energy physics, and be decoupled like 
other fermion doublers. This additional determinant may require more numerical cost than 
that of the gauge action, Eq. ()2.15|) . but it would be negligible compared to the cost of the 
dynamical overlap fermions. 
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3 Lattice simulations 



In this section, we explain our setups for the numerical simulations. 
3.1 Quenched QCD with admissible gauge fields 

Although several types of the gauge action that generate the "admissible" gauge fields sat- 
isfying the bound Eq. ()1.2|) are proposed [33,34,38], we take the simplest choice: Eq. fl2.15|) . 
We use three values of 1/e: 1, 2/3, and 0. Note that 1/e = corresponds to the conventional 
Wilson plaquette gauge action. The value 1/e = 2/3 is the boundary, below this value the 
gauge links can take any value in the gauge group SU{3) (all configurations are admissible.) 
and the positivity is guaranteed [93]. 

The link fields are generated with the standard hybrid Monte Carlo algorithm [25] (See 
appendix EI)- The molecular dynamics step size Ar is taken in the range 0.01-0.02 and the 
number of steps in an unit trajectory, Nmds, is 20-40. Every molecular dynamics step, we 
check whether the condition 1 — ReTrP^i,(x)/3 < e is satisfied or not. In fact, no violation of 
this bound was observed in our simulations. For thermalization, we discarded at least 2000 
trajectories before measuring observables. 

To generate topologically non-trivial gauge configurations, we use the initial link vari- 
ables, 



f/2(x) = I g2^ix35,2,t-i/L 

-27rix3<5^ i_i/L 




UJx) 



1 

^ g27rixiQinit5a;4,L-l/i' j 



(3.1) 



which is a discretized version of a classical solution with the topological charge Q = Qinit 
on a four- dimensional torus [94]. We confirmed that the topological charge assigned in this 
way agrees with the index of the overlap operator with s = 0.6. 

The simulation parameters and the plaquette expectation values (for the run with the 
initial configuration with Qinit = or 1) are summarized in Table ^ The length of unit 
trajectory is 0.07-0.4, and the step size is chosen such that the acceptance rate is larger than 
~ 70%. 
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Lattice 


1/e 


P 


At 


^mds 


acceptance 


Qinit 


plaquette 




1 


1.0 


0.01 


40 


89% 





539127f9) 


quenched 




1.2 


0.01 


40 


90% 





0.566429(6) 






1.3 


0.01 


40 


90% 





0.578405(6) 




2/3 


2.25 


0.01 


40 


93% 





0.55102(1) 






2.4 


0.01 


40 


93% 





0.56861(1) 






2.55 


0.01 


40 


93% 





0.58435(1) 







5.8 


0.02 


20 


69% 





0.56763(5) 






5.9 


0.02 


20 


69% 





58190f3) 






6.0 


0.02 


20 


68% 





0.59364(2) 


16^ 


1 


1.3 


0.01 


20 


82% 





0.57840(1) 


quenched 




1.42 


0.01 


20 


82% 





0.59167(1) 




2/3 


2.55 


0.01 


20 


88% 





0.58428(2) 






2.7 


0.01 


20 


87% 





0.59862(1) 







6.0 


0.01 


20 


89% 





59382f5) 






6.13 


0.01 


40 


88% 





0.60711(4) 


20^ 


1 


1.3 


0.01 


20 


72% 





0.57847(9) 


miPTiched 




1.42 


0.01 


20 


74% 





59165fl) 




2/3 


2.55 


0.01 


20 


82% 





0.58438(2) 






2.7 


0.01 


20 


82% 





0.59865(1) 







6.0 


0.015 


20 


53% 





0.59382(4) 






6.13 


0.01 


20 


83% 





0.60716(3) 




1 


0.75 


0.01 


15 


72% 





0.52260(2) 


with det 


2/3 


1.8 


0.01 


15 


87% 





0.52915(3) 







5.0 


0.01 


15 


88% 





0.55374(6) 


16^ 


1 


0.8 


0.007 


60 


79% 


+1 


0.53091(1) 


with det 


2/3 


1.75 


0.008 


50 


89% 


+1 


0.52227(3) 







5.2 


0.008 


50 


93% 


+1 


0.57577(3) 



Table 1: Simulation parameters and the plaquette expectation values. 
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3.2 Cooling method to measure the topological charge 

In order to measure the topological charge, we develop a new "cooling" method. It is 
achieved by the hybrid Monte Carlo steps with an exponentially increasing coupling, /3cooh 
and decreasing step size, Atcooi, as a function of trajectory, n^, i.e. 

/3cool = Anit X (l-S)"*, 

AtcooI = Arinit X (1.5)-"'/2^ (3.2) 

using the gauge action Eq. ()2.15|) with fixed 1/e = 1/ecooi- Note that a/AooiAtcooi is fixed in 
order to keep each step small. In this method, one can "cool" the configuration smoothly, 
keeping the admissibility bound, Eq. ()1.2|) . with 1/e = l/ccooi- The parameters are chosen 
as (Anit, Arinit, 1/ecooi) = (2.0, 0.01, 1) for the configurations generated with 1/e = 1, and 
(3.5, 0.01, 2/3) for the configurations with 1/e = 2/3 or 1/e = 0. Even for the gauge fields 
generated with the Wilson plaquette gauge action (1/e = 0), the condition 1/ecooi = 2/3 
can be used because it allows all values of SU(3). The link variables are cooled down after 
50-200 steps close to a classical solution in each topological sector. Then the geometrical 
definition of the topological charge [95], 

Qseo ^i^H e^^'^'^ReTr (P^,(x)P,.(x)) , (3.3) 

for the "cooled" configuration gives a number close to an integer times a universal factor 
Z~}^, namely, Q = ZgeoQgeo is close to an integer. The value Zgeo = 0.923(4) is calculated 
through would-be Q = I gauge configurations. As Fig. [T] shows, the topological charge Q is 
consistent with the index of the overlap-Dirac operator with s = 0.6, which is measured as 
described in Section 13.41 The consistency seems better for 1/e = 1 than for the standard 
Wilson plaquette action (1/e = 0). 

3.3 Large negative mass Wilson fermion 

Now let us discuss the lattice QCD with the determinant of Wilson fermion with negative 
mass at cut-off scale —(1 + s)/a as seen in Eq. fj2.18p . This is a more direct way to suppress 
the small eigenvalue of Hw, or keep topological charge, since the determinant det a^H^ 
would give a very small Boltzmann weight if Hw touches zero. For the gauge action, we use 
the topology conserving action Eq. (|2.15p again with three values of 1/e: 1, 2/3, and 0. We 
set the hopping parameter Hhop = 1/(8 — (1 + s)) = 0.208, which is known to be an optimal 
choice for the locality of the overlap Dirac operator [96], and the conventional pseudo- fermion 
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L=16, beta=1.42, 1/e=1 



L=16, beta=6.13, 1/e=0 





the overlap index — i — 




geometorical Q « 
















2110 2120 2130 2140 2150 2160 2170 2180 2190 2200 
conf 



8000 8100 8200 8300 8400 8500 8600 8700 8800 8900 9000 
conf 



Figure 1: Comparison of the topological charge Q; the geometrical definition Q = ZgeoQgeo 



after "cooling", and the index of the overlap-Dirac operator with s = 0.6. The agreement is 
better for 1/e = 1 (left) than for 1/e = (right). 

method is performed in the hybrid Monte Carlo steps to calculate the fermion determinant. 
For thermalization, we performed more than 500 trajectories before measuring observables. 

The topological charge is evaluated in the same way as explained in the previous sub- 
section, namely, by the geometrical definition Eq. (|3.3|) after 100-200 of the quenched HMC 
steps with increasing (3, switching off the Wilson fermion. Although the topological charge 
is expected to be more stable in this method, there may be a large scaling violation. Careful 
comparison with quenched studies without the determinant would be important. 

The simulation parameters are summarized in Table. Q 

3.4 Numerical implementation of the overlap Dirac operator 

For the overlap Dirac operator, Eq. ()2.1|) . we use the Wilson Dirac Hamiltonian aHw = 
'j^{aD]Y — 1 — s) with s = 0.6, which is empirically known as an optimal choice to suppress 
the small eigen-modes of laif^/l at 1/a ~2.0GeV in quenched QCD. The sign function is 
approximated by the Chebyshev polynomial (See Fig. |21); 

sgncheb(«^VK) = aHw/sqTtck^^{a^H^), 

where the function sqrtQ^^^^a^H^) in the range [a,b] is written 

b + a 2 



(3.4) 



1 / sqrt 



Chcb 



-^poly 

1=0 



X 



-X 



V2 



COS ri' 



\J {h + a) — {h — a) cos 
t = cosd). 



cosn- 



(3.5) 
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Note that the Chebyshev basis Tj(x) (of order i) satisfies the orthogonal relation 



1 



dt 



:T,(t)T,(t)=5,,. 



(3.6) 



1 



We use the numerical package ARPACK [97], which implements the implicitly restarted 
Arnoldi method, to measure the lowest eigenvalue a = A^j^ and the highest b = X^.^^. In 
some cases, we calculate the 10 lowest modes explicitly, project them out of Hw, and use 
the approximation Eq. ()3.4j) after the projection, in the range [an, 6] where we use the 11th 
eigenvalue for the lower bound; an = Xf^. 

In order to determine the order of the polynomial, A^poiy, it is important to note that the 
accuracy is expressed as a function of A^poiy and the condition number k, = Amax/ Amin = yb/a, 




for a random vector \v), where the constants A ~ 0.3 and B ~ 4.2 are turned out to be 
independent of /5, e and the lattice size L, as seen in Fig. 01 Note that the Ihs. of Eq. (j3.7|) 
is zero when the sign function is exact and the approximation breaks down when 1/k = 0. 

We also use ARPACK to calculate the eigenvalues of P DP and P^DP^, where P± = 
(1 ± 75)/2. The index can be obtained from the number of zero-modes of chirally projected 
Dirac operators. To construct the non-zero mode, |A), of D, the formula 




{v\v) 



Aexp{-BNpo\y/K) 



(3.7) 



D-X* 



D-X* 



P_|A) 



(3.8) 



ImA 



ImA 



is useful. 
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Chebyshev approximation 




_-| 5 I I I I I I I I I 

-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 

X 

Figure 2: The Chebyshev approximation of sgn function. We set A^poiy = 60 and the range 
[a = 0.1,6= 1]. 



(1+s)=1.6, L=16 
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(1+s)=1.6, L=20 



0.01 
1e-04 
1e-06 
1e-08 
1e-10 
1e-12 
1e-14 
1e-16 
1e-18 





beta=1 .42,1/6=1 »< 




beta=1 .3,1/6=1 □ 




beta=2.7,1/e=2/3 + 




beta=2.55, 1/6=2/3 




beta=6.13, 1/6=0 v 




b6ta=6.0, 1/6=0 




0.3*6xp(-4.2*x) - 







4 5 6 7 

NpOL*HWio/Hw„a, 



Figure 3: The accuracy Eq. ()3.7|l as a function of N^oi/k on a 16^ lattice (left) and a 20^ 
lattice. We use 4 gauge configurations and 10 values of A^poi = 60-195 for each parameter 
set. 
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4 Wilson loops and the static quark potential 



In order to confirm the practical feasibility of the topology conserving actions, Eg . (12.151) 
and Eq. ()2.18p . a careful comparison with the Wilson plaquette action should be done. In 
this section, we study gluonic quantities or Wilson loops to determine the lattice spacing, to 
estimate the scaling violations due to the change in the actions, and to test the perturbation 
theory with tadpole improved coupling. Then we can judge whether our naive expectation 
that 1/e effect is small (of order 0{a^)) or the effect of the determinant, det H^r, is negligible, 
is true or not. 

4.1 The static quark potential 

In the following, we assume that the topology of the gauge field does not affect the Wilson 
loops when the lattice size is large, say L > 1 fm, and choose the run with Q = or Q = 1 
initial configuration for the measurement. 

Wilson loops, W{r,t)'s, are measured using the smearing technique according to [98], 
where the spatial separation f/a is taken to be an integer multiples of vectors v = (1, 0, 0), 
(1, 1, 0), (2, 1, 0), (1, 1, 1), (2, 1, 1) and (2, 2, 1). Assuming the Wilson loop is an exponential 
function 

{W{f, t)) = exp{-V{f)t), (4.1) 

for large t/a, we extract the static quark potential aV{f). The measurements are done every 
20 trajectories and the errors are estimated by the jackknife method. 

As a reference scale to determine the lattice spacing, we evaluate the Sommer scales tq 
and Tc [99,100] defined by rQF(ro) = 1.65 and r'^F{rc) = 0.65, respectively. The force F{r) 
on the lattice is given by a differentiation of the potential in the direction of u/a = (1, 0, 0); 

\u/a\ 

for r/|r| = (1, 0, 0). r/ is introduced to cancel 0{a) errors in the short distances at tree level, 

1 aG{f) — aG{f— u) 

A-K^ri/aY \u/a\ 



«G'(r) = / '^-l :L'' (4.3) 



d^k n^=i cos{rjkj/ a) 
.(27r)3 4E?=isin2(A;,/2) 

In Table 12] we present the values of the Sommer scales ro/a, r^/a, their ratio Tc/tq, and 
the lattice spacing a (We assume tq ~ 0.5fm.). The data of aV{r) and rjF{rj) in the case 
that f/a is an integer multiples of u/a are given in Tables 121 ID and El 
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Lattice size 


1/e 




statistics 


Tn/a 


Tela 




a 




1 


1.0 


3800 


3.257(30) 


1.7081(50) 


0.5244(52) 


0.15fm 


quenched 




1.2 


3800 


4.555(73) 


2.319(10) 


0.5091(81) 


0.1 Ifm 






1.3 


3800 


5.140(50) 


2.710(14) 


0.5272(53) 


O.lfm 




2/3 

/ 


2.25 


3800 


3.498(24) 


1.8304(60) 


0.5233(41) 


0.14fm 






2.4 


3800 


4.386(53) 


2.254(10) 


0.5141(61) 


0.1 Ifm 






2.55 


3800 


5.433(72) 


2.809(18) 


0.5170(67) 


0.09fm 


16^ 


1 


1.3 


2300 


5.240(96) 


2.686(13) 


0.5126(98) 


O.lfm 


quenched 




1.42 


2247 


6.240(89) 


3.270(26) 


0.5241(83) 


0.08fm 




2/3 


2.55 


1950 


5.290(69) 


2.738(15) 


0.5174(72) 


0.09fm 






2.7 


2150 


6.559(76) 


3.382(22) 


0.5156(65) 


0.08fm 


144 


1 


0.75 


162 


4.24(15) 


2.240(37) 


0.528(24) 


0.12fm 


with det if^ 


2/3 


1.8 


261 


4.94(19) 


2.361(26) 


0.478(19) 


O.lOfm 







5.0 


162 


4.904(90) 


2.691(42) 


0.549(13) 


O.lOfm 


16^ 


1 


0.8 


207 


4.81(17) 


2.442(48) 


0.508(20) 


O.lOfm 


with det 


2/3 


1.75 


189 


4.71(19) 


2.279(48) 


0.484(22) 


O.llfm 







5.2 


225 


7.09(17) 


3.462(55) 


0.489(13) 


0.07fm 



continuum hmit [100] 0.5133(24) 



Table 2: Sommer scales ro/a, rg/a and their ratio. Here we assume Tq ~0.5fm. 
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1/e = 1 
(5 r/a Ti/a 


12^ oupnclipd 
ay(f) rjF{ri) 


16^ mipnchpd 
aV{r} rjF{ri) 


1.0 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 


0.50459(20) 

0.77828(61) 0.5056(10) 
0.9629(15) 0.9520(69) 
1.1176(27) 1.691(26) 
1.2623(45) 2.751(80) 
1.4052(77) 4.33(22) 




1.2 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 


n 44877fl6^ 

0.65982(39) 0.38993(65) 
0.78291(80) 0.6346(34) 
0.8775(13) 1.034(10) 
0.9588(29) 1.545(45) 
1.0322(47) 2.23(12) 




1.3 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 

7 6.41 

8 7.43 


0.42730(10) 

0.61711(34) 0.35252(99) 
72140f69i 53909f48i 
0.7977(12) 0.848(14) 
0.8608(21) 1.240(36) 
0.9230(25) 1.887(85) 


0.42709(20) 

0.61710(66) 0.35099(68) 
72130f92i 5490f29) 
0.7961(15) 0.8325(81) 
0.8583(23) 1.180(32) 
0.9150(27) 1.809(79) 
0.9636(51) 1.93(24) 
1.0215(51) 3.09(37) 


1.42 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 

7 6.41 

8 7.43 




0.40443(15) 

0.57416(43) 0.31444(58) 
0.66091(75) 0.4567(22) 
0.7200(12) 0.6583(61) 
0.7691(17) 0.940(14) 
0.8076(24) 1.189(48) 
0.8457(30) 1.675(64) 
0.8832(37) 1.91(14) 



Table 3: Potential and force values in the case that f/ a is an integer multiples of the unit 
vector u/a — (1,0,0). Results with 1/e = 1 (quenched). 



19 



1/e = 2/3 
(5 r/a Ti/a 


12^* quenched 
aV{r) r]F{ri) 


16'* quenched 
aV{f) r'jFijri) 


2.25 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 


0.48470(15) 

0.74012(57) 0.47189(97) 
0.9077(13) 0.8640(56) 
1.0463(22) 1.515(21) 
1.1701(38) 2.353(64) 
1.2901(58) 3.64(15) 




2.4 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 


0.44908(12) 

0.66434(41) 0.39770(70) 
0.79152(84) 0.6557(37) 
0.8889(15) 1.065(12) 
0.9749(23) 1.635(32) 
1.0541(30) 2.401(74) 




2.55 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 

7 6.41 

8 7.43 


0.42013(11) 

0.60682(36) 0.34493(58) 
0.70826(72) 0.5230(28) 
0.7806(13) 0.7913(90) 
0.8430(18) 1.187(18) 
0.8986(23) 1.686(37) 


0.42042(16) 

0.60786(51) 0.34590(72) 
0.71227(95) 0.5337(32) 
0.7878(16) 0.8211(93) 
0.8538(22) 1.210(21) 
0.9157(29) 1.765(47) 
0.9710(43) 2.229(84) 
1.0266(52) 2.94(15) 


2.7 1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 

7 6.41 

8 7.43 




0.39590(15) 

0.56100(44) 0.30650(53) 
0.64733(62) 0.4456(22) 
0.70527(90) 0.6329(56) 
0.7528(14) 0.907(14) 
0.7937(19) 1.309(28) 
0.8321(24) 1.531(44) 
0.8703(29) 2.035(80) 



Table 4: Same as Table El but with 1/e = 2/3 (quenched). 
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1/6=1 

r/a Tj/a 


14^ detK^ /3 = 0.75 
aV{r) r'jF{rj) 


16^ det H'tv (3 = 0.8 
aV{r) r'jF{rj) 


1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 

7 6.41 

8 7.43 


0.50199(51) 

0.7198(18) 0.4025(31) 
0.8480(32) 0.661(14) 
0.9457(50) 1.068(49) 
1.037(55) 1.74(12) 
1.097(85) 1.83(26) 


0.48708(37) 

0.6943(13) 0.3829(22) 
0.8137(25) 0.616(14) 
0.8936(42) 0.873(45) 
0.9681(57) 1.416(87) 
1.0334(85) 1.98(28) 
1.110 (11) 3.01(37) 
1.151(14) 2.26(80) 




1/6 = 2/3 

r/a rj/a 


14^ detH'^ (3 = 1.8 
aV{f) rjF{ri) 


16* detH'^ (3 = 1.75 
aV{r) r]F{rj) 


1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 

7 6.41 

8 7.43 


0.49349(42) 

0.7125(15) 0.4046(25) 
0.8323(25) 0.6178(91) 
0.9294(48) 1.062(40) 
1.0000(68) 1.34(12) 
1.0655(92) 1.99(24) 


0.50630(46) 

0.7363(18) 0.4249(31) 
0.8605(35) 0.640(20) 
0.9680(68) 1.175(65) 
1.0423(88) 1.41(16) 
1.118(15) 2.30(39) 
1.185(17) 2.64(55) 
1.261(24) 4.1(1.2) 




1/6 = 
r/a Tj/a 


14^ det H'^ (3 = 5.0 
aV{r) r]F{rj) 


16* det H'^ (3 = 5.2 
aV{r) rjF{ri) 


1 

2 1.36 

3 2.28 

4 3.31 

5 4.36 

6 5.39 

7 6.41 

8 7.43 


0.45795(42) 

0.6547(12) 0.3634(19) 
0.7606(22) 0.5466(84) 
0.8369(32) 0.833(28) 
0.9000(51) 1.200(70) 
0.9730(75) 2.21(13) 


0.42080(26) 

0.58846(72) 0.3098(13) 
0.6722(12) 0.4318(48) 
0.7246(16) 0.617(14) 
0.7745(25) 0.872(43) 
0.8120(27) 1.137(60) 
0.8426(37) 1.21(10) 
0.8762(39) 1.84(15) 



Table 5: Same as Table 01 with det 
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Tc/tq is a good quantity to estimate the scaling violation, comparing with the value 
^c/'^o = 0.5133(24) in the continuum limit, which was obtained with the plaquette action 
[100]. Fig. 13 shows the dependence of this ratio with different values of 1/e. The quenched 
results with 1/e = 2/3 and 1 agree very well with Ref [100] except for the coarsest lattice 
points around (a/rg)^ ~ 0.1. Also in the case with detif^, as seen in Fig. El any large 
scaling violation has not been seen although the statistics are much poorer. 

Fig. ini presents the quark potential itself in a dimensionless combination, i.e. V{r) = 
f'o{y{'f^) — yi'f'c)) as a function of |r|/ro. V{rc) is evaluated by an interpolation of the data 
in the direction r/\r\ = (1,0,0). The data at /5 =1.3, 1/e = 1 in the quenched case and 
at /3 = 0.75, 1/e = 1 with detif^ are plotted together with the curve representing the 
continuum limit obtained in [100]. The agreement is satisfactory (less than two sigma) for 
long distances \r\/rQ > 0.5. For short distances, on the other hand, one can see deviations 
of order 10%, as shown in Fig. [3 where the ratio {V{f) — Kont(|^|))/Kont(|^) is plotted. 
V'contd'l) represents the curve in the continuum limit drawn in Fig. El The plots at r/a = 
(1, 0, 0) and (2, 0, 0) obviously deviate from zero in the upward direction, while the plots 
at (1,1,0) and (1,1,1) are lower than zero. This fact indicates the rotational symmetry 
violation in a short distance. Let us concentrate on the point (1, 0, 0) as a function of the 
lattice spacing, as seen in Fig. |S1 (left panel). We find that the size of the violation is quite 
similar to that with Wilson plaquette action and independent of 1/e. One might think that 
the rotational symmetry goes worse in the continuum limit, but one should note the fact 
that the relevant scale of the observable also diverges as 1/a. After a correction at the tree 
level by introducing dj as 1/(4716?/) = G{d), which is an analogue of rj in Eq. (j4.3p but is 
defined for the potential, we obtain the plot on the right panel of Fig. |H1 where one can see 
an improvement. The remaining deviation should be of order ^^(l/a), which would vanish 
as ~ 1/ ln(l/a) near the continuum limit. The data at r = (1, 0, 0) with det H^r have larger 
statistical errors but show a very similar behavior. 

These observations show that both of scaling violations and rotational symmetry vi- 
olations are reasonable and the topology conserving actions are feasible in the numerical 
studies. 

Finally, we confirm our assumption that the topology does not affect the quark potential 
by measuring tq for two initial value of Q (0 and —3). Measurements are done on a 16^ 
lattice at (3 = 1.42, 1/e = 1, for which the probability of the topology change is extremely 
small as discussed in Sec. El Our results are rg/a = 6.24(9) for the Q = initial condition 
and 6.11(13) for Q = —3. Details are seen in Table. El 
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0.49 
0.48 
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0.06 



0.08 



0.1 



Figure 4: Quenched results (without det H^) oirdrQ. Squares and triangles are data for the 
topology conserving gauge action with 1/e = 2/3 and 1, respectively. Open circles represent 
the Wilson plaquette gauge action (from [100]) and the filled circle is their continuum limit. 
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Figure 5: The same as FigEJwith det H\ 



w 



Table 6: The quark potential with different initial topological charge. The results are at 
1/e = 1.0 and (3 = 1.42 on the L = 16 lattice. Any large discrepancies are not seen. 



Qinit Ar 


^mds 


acceptance 


Stabg 


plaquette 


ro/a 




0.01 


20 


82% 


961 


0.59167(1) 


6.240(89) 


0.5241(83) 


-3 0.01 


20 


83% 


514 


0.59162(1) 


6.11(13) 


0.513(12) 
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beta=1.3 1/e=1 beta=0.75, 1/e=1 with det H^,,^ 




Figure 6: Static quark potential at /3 = 1.3, 1/e = 1 (quenched) on a 12'^ lattice (left) 
and P = 0.75, 1/e = 1 (with detif^) on a 14^ lattice (right) . The curve is the continuum 
limit obtained by an interpolation of the data of [100]. Different symbols show V^(r)'s with 
different orientations parallel to u's. 
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Figure 7: Violation of rotational symmetry (y{r) — VJ;ont(|?1))/Kont(|?1), where Kontd^) 
denotes the continuum limit. Results for P = 1.3, 1/e = 1 (quenched) and P = 0.75, 1/e = 
1 (with det Hy^) are shown. The error of Kontd^) is not taken into account (< 1%). 
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Figure 8: (V"(r) — Kont(|^))/Kont(|^) for r = (1,0,0). Quenched results with different /? 
and 1/e values are plotted as a function oi a/vQ. Left and right panels show the plot without 
and with the tree level improvement of the argument r, respectively. See the text for details. 



4.2 Perturbative renormalization of the gauge coupling 

From 1x1 Wilson loop, W{1, 1), the so-called "mean field improved" bare gauge coupling 
can be defined. In quenched study with the standard Wilson plaquette action, it is known 
that the perturbation theory with mean field improvement coupling converges very well. 
Here we would like to define the mean field improved coupling with the topology conserving 
gauge action, Eq. ()2.15|) . evaluate the coupling renormalization and see its convergence in 
2-loop perturbation theory. 

Two-loop corrections to the gauge coupling for general one-plaquette actions (constructed 
by the plaquette only and no rectangular term is involved.) is calculated by Ellis et al. [101]. 
Using their formula, the renormalized gauge couping qm defined in the "Manton" scheme is 
given by 



9 m r 

where g denotes the bare coupling and the coefficients Ai and A2 are calculated as 



(4.4) 



A, 



A, 



2M - 3 



S4- 



+ U{Nl + 1) 



a^[s,{2Nl - 3) + UN,{Nl + 1)] + s,'-^^^^^l^^f±^ 
— + h-{Nl + l){Nl + 3) 



—s 



8N, 

2iv! 



— 254^4 



(2iV^^-3)(iV^^ + 2) 

Nr. 



-t^(iVe' + l)(A^c+2). (4.5) 
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Here, the parameters are A^'^ = 3, S4 = -1/4!, sq = 1/6!, U = l/(4A/'ce), te = l/iSN^e"^), 
uq = — l/(4!Ai"ce), and = —0.0011(2). The values of the next-to-leading and next-to-next- 
to-leading order coefficients Ai and A2, when 1/e = 1,2/3 and 0, are given in Tabled 
Let us define the mean field improved coupling by 

I + 'l-^V^ ] (46) 

with the measured value of the plaquette expectation value P = {W{1, 1)) (see Table [T]). It 
is defined from a coefficient of when we rewrite Pf^u = -Pexp(ia^F^i,) and expand the 
action Eq.^J^ [102] . 

A perturbative expectation value of the plaquette, or 1 x 1 Wilson loop, is evaluated with 
the general one-plaquette action by Heller et al. [103] as 

{W{l^)) = 1-/ ^ % ,(l,l)-^^(Ar^^-l)X(l,l) 

Here W^2(l, 1) and X(l, 1) are from the original calculation [104] for the Wilson plaquette 
gauge action, and Z(l, 1) = (1 — 1/1^)14^2(1, l)/4 (on a symmetric lattice V = L^) is in- 
troduced for generalization. Their values are ^2(1,1) = l/^j -^(1? 1) = —1-01 x 10~^ and 
Z{1,1) = 1/32 in the infinite volume limit. The constant C is written 



C 



(4. 



where C2{R) denotes the quadratic Casimir operator in a representation R of the group 
SU{Nc). dji is the dimension of the representation R, and T{R) is defined such that 
TinitH^) = T{R)6°'^ for the group generator t". The coupling sr{P) is defined as a cou- 
pling when we rewrite the gauge action in terms of a general form of the one-plaquette 
action, 

= E E 'Ri(^) fl - ;^ReTr^P« (x)j , (4.9) 

where denotes the plaquette variable in the R representation. The values of these 
parameters for the topology conserving gauge action Eq. ()2.15p are 

S3(/?) = (1 + /?, se{f3) = -^^(3, s8(/?) = -^/?, (4.10) 

The other parameters are T(3) = 1/2, T(6) = 5/2, T(8) = 3, C2(3) = 4/3, C2(6) = 10/3 
and 6*2(8) = 3. With these numbers, we obtain C = 5 — 20/e and we finally get 
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A2 


Bi 


B2 




A2-B2 





-0.20833 


-0.03056 


-0.33333 


-0.03472 


0.12500 


0.00416 


2/3 


0.34722 


-0.04783 


0.11111 


-0.05015 


0.23611 


0.00233 


1 


0.62500 


-0.10276 


0.33333 


-0.13194 


0.29167 


0.02919 



Table 7: Next-to-leading and next-to-next-to-leading order coefficients of the coupling renor- 
malization for various e. Mean field improved coefficients Ai — Bi, A2 — B2 are also presented. 
See the text for details. 

Since the perturbative relation between the bare coupling g and the boosted coupling g 
Eq. (|4.6|) is given by 

4 = ^ + ^1 + ^2/, (4.12) 
9^ 9 

where 

3V e/' V e/ V18e 144/ 9e 3e2' ^ ^ 

one obtains the Manton coupling in terms of the boosted coupling as 

^ = 4 + (^1 - 5i) + {A2 - B2)f. (4.14) 
9m 9^ 

Numerical values of Bi and Ai — Bi are listed in Table [7| One can see the effect of the mean 
field improvement; the two- loop coefficient A2 is significantly reduced by reorganizing the 
perturbative expansion as in Eq. ()4.14|) . 

With these results, the Manton coupling is obtained for each lattice parameters. In 
Fig. ini we plot the coupling evaluated at a reference scale 5/ro as a function of a? /r"^. We 
use the two- loop renormalization equation to run gM to the reference scale. Although gM 
is very different at the tree level, the one-loop results are already in good agreement among 
the different values of 1/e. The two-loop corrections show that the perturbative expansion 
converges very well and a good agreement among different 1/e. Moreover, the scaling toward 
the continuum limit seems also good in the two-loop level plots. 

5 Stability of the topological charge 

5.1 Admissibility condition and topology stability 

As discussed in Sec|21 it is difficult to construct an exact and practical geometrical definition 
of the topological charge for non-Abelian theories, unless the gauge fields are very smooth 
[105] (note that ()3.3|) gives non-integers). But even if we choose larger e ^ 1/20, it is quite 
possible that the barriers among the topological sectors are high enough to suppress topology 
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Figure 9: l/f?!/ in the Manton scheme are plotted. The mean field improved expression 
Eq. (|4.14| ) is used with the measured plaquette expectation value. Different symbols distin- 
guish the value of 1/e (1 for circles, 2/3 for squares, for triangles). Open symbols with 
dot-dashed line represent tree-level results and open symbols with dashed lines are one-loop. 
The best results including two-loop corrections are shown by filled symbols with solid lines. 



changes for hundreds of HMC trajectories, since any gauge action has a tendency to prevent 
the topology transitions in the continuum limit. 

In Table. |H1 we present our data of the topological charge stability. 



Stab, 



Q 



trj 



1"plaq X #(5 ' 



(5.1) 



where Xpiaq is the plaquette autocorrelation time, which is measured according to Appendix 
E of [106]. A'^trj denotes the total length of the HMC trajectories and j^Q is the number 
of topology changes during the trajectories. This definition, Stabg, gives a mean number 
of uncorrelated gauge configurations which can be sampled without changing the topology 
along the simulations. But we should mention that it only gives an upper limit, because the 
topological charge is measured only once per 10-20 trajectories and we may underestimate 
the topology changes if Q jumps and returns quickly to the original value within this interval. 
Therefore, Stabg is not reliable when the topology change is very frequent. 

In Fig.Hnithe results are plotted as a function of the lattice spacing squared. Clearly, one 
can see that the stability goes better as 1/e increases when the lattice spacing is the same. 
Also, the stability gets worse as the lattice size is increased from L/a = 12 to L/a = 20. This 
is no surprising, because the topology change occurs through local dislocations of gauge field 
and its probability scales as the volume. An important notice here is that the topological 
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Lattice 


1/^ 




To/ a 




^plaq 


#Q 
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1.0 


3.257(30) 
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2.91(33) 
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1 o 
i .z 


'±.000[ (6) 


1 snnn 
ioUUU 


L.Oi)[LO ) 


zoo 


4o 






1.3 


5.140(50) 


18000 


1.091(70) 


69 
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2/3 


2.25 


3.498(24) 
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5.35(79) 
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2.4 


4.386(53) 
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2.62(23) 
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2.55 


5.433(72) 
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2.86(33) 


123 
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1.3 


5.240(96) 
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3.2(6) 


78 


46 
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1.42 


6.240(89) 


5000 


2.6(4) 


2 


961 




2/3 


2.55 


5.290(69) 


12000 


6.4(5) 
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18 






2.7 
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14000 


3.1(3) 


6 


752 
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14 
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1.42 
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7000 


3.8(8) 


29 
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Z.OO 
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2.7 





7800 


3.5(6) 


20 
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1600 


14.4(7.8) 
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3 






6.13 
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1298 


9.3(2.8) 


4 


35 


14^ 


1 


0.75 


4.24(15) 


3500 


5.05(82) 





>693 


with det 


2/3 


1.8 


4.94(19) 


5370 


11.1(2.1) 





>483 







5.0 


4.904(90) 


3120 


21.4(6.5) 





>146 


16^ 


1 


0.8 


4.81(17) 


627 


0.69(10) 





>908 


with det 


2/3 


1.75 


4.71(19) 


580 


1.64(40) 





>353 







5.2 


7.09(17) 


730 


1.54(27) 





>474 



Table 8: The stabihty of the topological charge, Stabg. The HMC trajectory length A^trj, 
the plaquette autocorrelation time Tpiaq, and the number of topology change -^Q are also 
summarized. 
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Figure 10: Stability of the topological charge in the quenched study. The results of different 
lattice size are plotted with different symbols; L/a = 12 (circles), 16 (squares), 20 (upward 
triangles). The value of 1/e is distinguished by the line type: 1/e = (dotted), 2/3 (dashed), 
1 (solid). Downward triangles are the data of [34] measured on a 16^ lattice. 



charge can be very stable in any case in the continuum limit. The stability quickly rises as 
a ^ 0. 

To study QCD in the e-regime in a fixed topological sector, the lattices (1/e,/?, L) ~ 
(1,1.42,16) and (2/3,2.7,16) would be appropriate. Their physical size is around L ~ 
1.25 fm and the topological charge is stable for (100 — 1000)rpiaq trajectories. 



5.2 Negative mass Wilson fermion to fix topology 

In contrast to the case with topology conserving gauge action Eq. (j2.15|) . one expects that 
the fermion determinant, detif^, can rigorously fix the topology of gauge fields along the 
simulation. In fact, as Table. |H1 shows, the topology change has never occurred in every run 
with different parameters, which are chosen such that the lattice spacing is around a ~ 0.08- 
O.lfm. Our data show that Q is unchanged for, at least, 100-1000 Tpiaq trajectories, even if 
one ignores the thermalization steps. 
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6 The effects on the overlap Dirac operator 



6.1 Low- lying mode distribution of Hw 

As explained in Sec 13 .41 the order of Chebyshev polynomial approximation, A^poiy has to 
be proportional to the condition number, k = Amax/ Amin, in order to keep a certain desired 
accuracy. Since both topology conserving actions, Eq. ()2.15|) and Eg. ()2.18j) would play a role 
to suppress the occurrence of low-lying eigenvalues of |aifiy| = \ j5{aDw — 1 — they may 
be useful to reduce the numerical cost of calculating the overlap Dirac operator. 

Fig. ^2 shows a typical comparison of the eigenvalue distribution on a 16^ lattice. The 
values of (3 and e are chosen such that the Sommer scale r^/ais roughly equal to 5 (left panel) 
and 6.5-7 (right panel) , which correspond to a ~ 0.1 fm and a ~ 0.08 fm respectively. From 
the plot we observe that the density of the low-lying modes is relatively small for larger 
values of 1/e in quenched case. Let us discuss this more quantitatively. In Table El we list 
the probability, P\^i^^<o.i, that a configuration has the lowest eigenvalue lower than 0.1. For 
the above example (left panel), the probability is 74% for the standard Wilson gauge action 
(/5 = 6.0, 1/e = quenched ), but it decreases to 53% (47%) for 1/e = 2/3 (1) (quenched). 
It is interesting to note that the data with det show that the average of Amin is not so 
large, but the occurrence of very small eigenvalues <^0.1 is strongly suppressed. The lowest 
eigenvalue in these configurations is 1.2e-03 in the quenched case at /9 = 6.0 and 1/e = 0, 
while it is 0.042 in the case with det if^ at /5 = 0.8 and 1/e = 1. 

Also, for another lattice spacing (ro/a — 6.5) and lattice size 20^, a similar trend can 
be seen. In Table M we also present the ensemble average of the lowest eigenvalue Amin and 
the inverse of condition numbers Amax/Amin and Amax/Aio, where Aio and Amax denote the 
10th and the highest eigenvalues respectively. We may conclude that the lowest eigenvalue is 
higher in average when we use the topology conserving actions. In the numerical implemen- 
tation of the overlap-Dirac operator, the low-lying eigen-modes of aH^ are often subtracted 
and treated exactly. Then the higher mode contributions are approximated by some poly- 
nomial or rational functions. Here, we assume that 10 lowest eigen-modes are subtracted 
and compare the relative numerical cost on the gauge configurations with different values of 
1/e. From Table El we observe that the reduced condition number is about a factor 1.2-1.4 
smaller for 1/e = 1 than that for the standard Wilson gauge action. 

We also check that the above observation does not change when one shifts the value of s 
in a reasonable range. In Fig. ^1 a typical distributions of the low-lying eigen-modes for s 
= 0.2-0.7 are plotted. We find that the advantage of the topology conserving actions does 
not change. Also, from these plots we can see that s ~ 0.6 is nearly optimal for all the cases. 



31 



10 lowest eigenvalues of H^j (L=16) 



10 lowest eigenvalues of H^,, (L=16) 
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Figure 11: Ten lowest eigenvalues of |a-ffvK| = |75(a.-Dvy — 1.6) | for gauge configurations with 
ro/a ~ 5 (left) and ro/a — 6.5-7 (right). Quenched data are shown for 1/e = 1 (pluses), 2/3 
(stars), (squares) and triangles are the result with detif^. The lattice size is 16^. 



lattice size 


lA 


P 


ro/a 




^min 


-^min/-^max 


Aio/ Amax 


20^ 


1 


1.3 


5.240(96) 


0.64 


0.0882(84) 


0.0148(14) 


0.03970(29) 


quenched 


2/3 


2.55 


5.290(69) 


0.75 


0.0604(53) 


0.0101(08) 


0.03651(27) 







6.0 


[5.368(22)] 


0.97 


0.0315(57) 


0.0059(34) 


0.02766(46) 




1 


1.42 


6.240(89) 


0.22 


0.168(13) 


0.0282(21) 


0.04765(32) 




2/3 


2.7 


6.559(76) 


0.19 


0.151(11) 


0.0251(19) 


0.04646(37) 







6.13 


[6.642(-)] 


0.45 


0.0861(83) 


0.0126(15) 


0.03775(50) 


16^ 


1 


1.3 


5.240(96) 


0.47 


0.111(12) 


0.0187(21) 


0.04455(31) 


quenched 


2/3 


2.55 


5.290(69) 


0.53 


0.1038(98) 


0.0174(16) 


0.04239(36) 







6.0 


[5.368(22)] 


0.74 


0.0692(90) 


0.0116(15) 


0.03451(62) 




1 


1.42 


6.240(89) 


0.07 


0.219(13) 


0.0367(21) 


0.05233(26) 




2/3 


2.7 


6.559(76) 


0.13 


0.191(12) 


0.0320(19) 


0.05117(29) 







6.13 


[6.642(-)] 


0.27 


0.139(10) 


0.0232(17) 


0.04384(38) 


16^ 


1 


0.8 


4.81(17) 


0.12 


0.1502(87) 


0.0255(11) 




with det 


2/3 


1.75 


4.71(19) 


0.55 


0.0999(44) 


0.0170(10) 









5.2 


7.09(17) 


0.05 


0.1771(69) 


0.0297(12) 





Table 9: The probability, -PAmin<o.i) to sample a configuration such that Amin < 0.1 where Amin 
is the lowest eigenvalue of the hermitian Wilson-Dirac operator [a-ffv^l = I'JbiO'Dw — 1-6) |. 
Averages of the lowest eigenvalue and the inverse of condition numbers are also listed. The 
Sommer scale ro/a is the results from L = 16 lattice simulations. The values with [■ ■ •] are 
from [100] with an interpolation in /3. 
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Figure 12: Ten near-zero eigenvalues on 16^ lattices with r^/a ~ 5.3 (left) and ro/a ~ 6.5 
(right). Results are plotted as a function of 1 + s. Plots for 1/e = 2/3 and are slightly 
shifted for clarity. 



6.2 Locality 

When we talk about the locality of the overlap-Dirac operator, it means that the norm 
\\D{x^y)v{y)\ \ with a point source vector v at Xq should decay exponentially as a function 
of |x-a;o| [29] 

\\D{x,y)v{y)\\r^CeM-D\x-x^\), (6.1) 

where C and D are constants. Actually we observe this property, as seen in Fig. E| The 
plot shows the results for different values of 1/e at the lattice scales ro/a — 5.3 (left) and 
6.5 (right). We find no remarkable difference on the locality when we change 1/e. 

Recently, it has been indicated that the mobility edge is more crucial quantity which gov- 
erns the locality of the overlap-Dirac operator [69-73]. It would be interesting and important 
to see the dependence of the mobility edge on the parameters in the topology conserving 
actions, which should be done in the future works. 
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Figure 13: {\\D{x,y)v{y)\\/\\D{0,y)v{y)\\)'^ with xq = measured on 10 gauge configurations 
for different values of 1/e (quencfied). Tlie lattice scale is ro/a ~ 5.3 (left) and 6.5 (right). 

7 Lattice QCD in the e-regime with fixed Q 

In previous sections we have discussed the topology conserving actions which might be helpful 
to the full QCD simulation with the overlap fermion determinant. Here we would like to 
show great uses of the overlap Dirac operator itself to measure physical observables. Once 
a set of configurations of Nf = 2 or 3 full QCD is generated, it would be much easier and 
more reliable to extract physical quantities such as pion mass, decay constants, Kaon bag 
parameters, and so on, using the overlap Dirac operator which respects the "exact" chiral 
symmetry. One of attractive numerical applications would be the so-called e-regime, where 
both of chiral symmetry and topology have significant effects on the observables. 

In this section we demonstrate the analysis of quenched lattice QCD in the e-regime to 
determine the low energy constants of the quenched chiral perturbation theory (qChPT). Of 
course, both quenched QCD and qChPT are unphysical and the data of pion decay constant 
or the chiral condensate have little relation to the true value of them in nature. However, this 
quenched study is still interesting to see how these low-energy constants can be extracted 
when the full QCD simulation is established in the future. 



7.1 Meson correlators in the e-regime 

It is believed that the low energy limit of QCD is described by the pion effective theory, or 
(q)ChPT. To determine the fundamental parameters of (q)ChPT is one of relevant issues 
of lattice QCD. In the so-called e-regime [43-48], where the size of the space-time box is 
smaller than the pion Compton wave length; L <^ I/ttItt (but larger than the QCD scale 
1/Aqcd L, which should be guaranteed so that the pion can be treated as a point 
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particle and other heavier hadrons are decoupled.), (q)ChPT with a expansion parameter 
~ itlt^/K ~ p^/A^, is still applicable. Here A is a cutoff scale of the chiral Lagrangian 
roughly around 1 GeV. An important notice is that the low energy constants in (q)ChPT 
are defined at the cutoff scale and it does not depend on whether the system is in the e- 
regime or the large volume regime. Therefore, once or S are determined in the e-regime, 
one can use them in the standard (q)ChPT in a larger volume. Analytic calculations of 
meson correlation functions in the e-regime have been widely studied for both ChPT and 
qChPT [107-109]. 

In the quenched case, the fundamental parameters of the effective theory are the pion 
decay constant, F^r, the chiral condensate, S, the singlet mass, mo, and the coefficient of an 
additional singlet kinetic term, a. Here we just present some results of 1-loop calculation, 
which are relevant to our numerical study. The details of calculation and the definition of the 
functions and coefficients are summarized in Appendix O The scalar condensate at 1-loop 
level in Q topological sector is given by 



- (tA^) = Sq(/x') = Eeff/i'(/|Q|(/iO^|Q|(/^0 + /|Q|+l(/i')^|Qhl(/iO) + Seff 



M 



(7.1) 



The triplet axial vector, scalar and pseudo-scalar correlators are given by 
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The singlet scalar and pseudo-scalar correlators are 
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A{x)b 
A{x)b+ 



(7.2) 
(7.3) 

(7.4) 



a+ + a_ + 4' 



(7.5) 
(7.6) 



As mentioned in AppendixEl the results are valid only when |Q| is small, which is a particular 
restriction of quenching the fermion determinant. 

As these equations of qChPT show, the meson correlators in the e-regime should be 
quite sensitive to the topological charge and the fermion mass. These prominent m and Q 
dependences are used to evaluate the low energy constants, or namely, we fit our lattice data 
with the above equations, in which F^, S, mo and a are treated as free parameters. The 
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parameter S always appears associated with the quark mass m as mE is renormahzation scale 
and scheme independent. The value of S in the following analysis should be understood as a 
bare quantity in the lattice regularization at a scale 1/a. To relate them with the conventional 
scheme such as the MS scheme requires perturbative or non-perturbative matching, which 
is beyond the scope of this work. 

7.2 Lattice observables with the exact chiral symmetry 

We generate gauge link variables at /? = 5.85 and 1/e = in the quenched approximation 
on a L^T = 10^ x 20 lattice. The spacial length L of the box is about 1.23 fm. We employ 
the massive Dirac operator for the valence quark, 

D^=[l~—^D + m, (7.7) 

where the bare quark mass are chosen to be very small; am =0.0016, 0.0032, 0.0048, 0.0064, 
0.008, which corresponds 2.6-13MeV. The number of configurations for each topological 
sector is given in Table ITTH We analyze the gauge configurations in \Q\ < 3 sectors. 

When the exact inversion of the overlap-Dirac operator is needed, we use the techniques 
described in Ref. [53]; for a given source vector r], we solve the equation 

D^ij = V, (7.8) 

by separating the left and right handed components as = P-ip + P+4' solving two 
equations 

= {P^DlD,^P_)-'P_Dl^, (7.9) 
P+ij = {P+D^P+)-\P+Ti-P+D^P.i;), (7.10) 

consecutively. (The above equations apply to positive Q cases and the same procedure 
applies with a replacement P+ ^ P_ to negative Q cases.) The conjugate gradient (CG) 
algorithm is used to invert the chirally projected matrices with the low-mode preconditioning 
in which 20 lowest eigen-modes are subtracted. 

As one of special features of finite volume regime, a remarkable low-mode dominance in 
the chiral limit is expected [110-112] . Using ARPACK [97] again, we calculate 200 |Q| 
lowest mode eigenvalues and their eigenfunctions of the overlap Dirac operator. The index 
Q is calculated at the same time. Note that these eigenvalues cover more than 15% of the 
circle in the complex space of the eigenvalues of aD as Fig. El shows. Then, the inverse of 
the overlap operator is decomposed as 

Mow 

g (^_^„^,),^^,„ ...(xWfa) + AC„-(x.,). (7,11) 
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where Ai's are eigenvalues of D and fi(x)'s are their eigenvectors. We set A^iow = 200 + 
In the following analysis, we often use the low-mode approximation of the expectation value 
(We will denote (■ ■ ■)iow instead of a simple bracket.) by inserting the low- mode part only, 

An^(a;,Z/)iow = 71 /oAx I (7-12) 

~i (1 — am 2]Xi + m 



1=1 



to the propagators. This would be valid only when it is proved that the higher mode's con- 
tributions are sufficiently small or can be canceled in the combination of different operators. 
An great advantage of using D^^{x,y)iow instead of exact D^{x,y) is that the inversion at 
any x and y is obtained without performing the CG algorithm, so that one can easily average 
the correlators over the space-time arguments, like 



'5:vr.V'(^,t)^r,^(o,o)\ (E^r,V'(f,t + to)V'r.^(5^o,to)) , (7.i3) 

^ Mow ^0,io \ ^ 'low 



where Fj denotes some gamma matrix. This so-called low-mode averaging dramatically 
reduces the fluctuation of the correlators as shown in Fig. which is also reported in 
Refs. [57,113]. 



Table 10: Number of configurations in each topological sector. 

\Q\ 12 3 

# of confs. 20 45 44 24 
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Figure 14: Lowest 202 eigenvalues of the overlap-Dirac operator aD at beta=5.85 on a 10^20 
lattice with topological charge Q = —2. The eigenvalues cover a vr/S arc of the circle. 
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Figure 15: The pseudo-scalar correlator at am = 0.008 and \Q\ = 1. Filled symbols denote 
the data with the low-mode averaging, where we use {L/2Y x T/2 source points, while open 
symbols are not averaged. 
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Figure 16: Triplet axial-vector (left) and pseudo-scalar (right) correlators at am = 0.008 
and \Q\ = 1. The low- mode-approximated correlator is compared with the corresponding 
exact one. The 20-|-|(5| low-modes are not sufficient to evaluate the axial correlators but 
they are enough for the pseudo-scalar correlators at long distances. 



7.3 Numerical results 

7.3.1 from the axial-vector correlator 

First let us consider the axial- vector current correlator ()7.2|1 . which is most sensitive to 
and not contaminated by the parameters itlq and a. Here we do not use the low- mode 
approximation because 200-1- |Q| modes are not sufficient to estimate the total propagator, 
as the left panel of Fig^J shows. 

A naive definition of the axial vector current is A^(a;) = ■?/5(a;)757^(r"/2)^/'(a;), constructed 
from the overlap fermion field ip{x). Note that it is not the conserved current associated 
with the lattice chiral symmetry, and (finite) renormalization is needed to relate it to the 
continuum axial- vector current. We follow the method applied in Refs. [114,115] to calculate 
the Za factor non-perturbatively. Namely, we calculate 

ai:.-(Vo^g(x,t)P°(0,0)) 

E.(P»(^,t)P»(o,o)) ' ^'-'^^ 

where Vo denotes a symmetric lattice derivative and P"(x) = ^/'(r"/2)75'?/' ^. 

Fitting the average of Pm{t) over all topological sectors with a constant in the range 

^We should use the chiraUy improved operator, P°^^p{x) — 'tp{T°- /2)j^{l — ^D)ijj, but for the on-shell 
matrix elements such as the one considered here, one can use the equation of motion to replace the aD by 
~am/{l — am/2), which is negligible for our quark masses. We therefore use the local operator for P°-{x). 
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Figure 17: ap as a function of the bare quark mass. Dashed hne is a result of the quadratic 
fit: ap = 0.00001 + 1.390(ma) - 0.19(ma)2. 



7 < t < 13, we obtain 



ap{ma) 



a(V^A^(x)P"(0)) 2ma 



(7.15) 



(P-(x)P-(O)) Za 

at four quark masses am = 0.0016, 0.0048, 0.008, 0.016, as plotted in Fig. El which shows 
a very good chiral behavior. With a quadratic fit we obtain 



ap{ma) = 0.00001(2) + 1.390(14)(ma) - 0.19(74)(ma) 



(7.16) 



The constant term is perfectly consistent with zero and we extract an accurate value, 
= 1.439(15) , which is consistent with the value Za = 1.448(4) reported in Ref. [116] which 
was done with the same (3 and s. 

Now one can compare the renormalized axial-vector correlation function with the qChPT 
result 



2Zl E(^o(x, t)Ao{0, 0)f = 2 ( ^ + 2mS|Q|(^)r/^i(|t/r|) ) . 



(7.17) 



to extract P^- and S. In Ref. [55] it is reported that the correlators suffer from large statistical 
fluctuation when |(5| = and the data in the other topological sectors are insensitive to E, 
but it turned out that two-parameter (F^ and S) fitting does work well when we treat the 
data of different topology and fermion masses simultaneously. As shown in Fig. our 
data at am = 0.0016, 0.0048, 0.008 in the \Q\ < 1 sectors are well described by the qChPT 
formula (j7.17|) . A simultaneous fit in the range 7 < t < 13 yields F^^ = 98.3(8.3) MeV and 
^1/3 _ 259(50) MeV with x^/dof = 0.19. The result for P^- is in a good agreement with that 
of the previous work [57], 102(4) MeV. 



40 



o 
II 

O 



< 

V 



O 




10 

t 



15 



20 



O 

CO 



< 

o 

< 



< 

V 



0.003 
0.0025 

0.002 
0.0015 

0.001 
0.0005 


0.003 
0.0025 

0.002 
0.0015 

0.001 
0.0005 




qChPT F=98MeV,I=(259MeV)'^ 



10 

t 



15 



qChPT F=98MeV,I=(259MeV)'^ 



10 
t 



15 



20 



20 



Figure 18: Axial-vector current correlators at am = 0.0016 (left) and 0.008 (right) for \Q\ 
= (top) and 1 (bottom). The dashed lines are the result of simultaneous fitting of the data 
for Q = and 1 at am = 0.0016, 0.0048, 0.008 in the region 7 < t < 13. 



On the other hand, the correlators at \Q\ = 2 do not agree with the above fit parameters 
as shown in Fig. [THl As discussed before, it may indicate that the topological sector \Q\ = 
2 is already too large to apply the qChPT in the e-regime. 

7.3.2 S, Seff and a from connected S and PS correlators 

We find that the scalar and pseudo-scalar triplet correlators are approximated precisely with 
the lowest 200-|-|Q| eigen-modes at small quark masses (m = 2.6-13 MeV) as the right 
panel of Fig^J shows. In the range 7 < t < 13, the systematic error of this low-mode 
approximation is estimated to be only 1% for the scalar correlators in \Q\ > 1 sectors and 
the pseudo-scalar correlators in all the topological sectors. In this way we measure 

(5(t))« = -2 X:(l + snS'ix + Xo)S'{xo))t^, (7.18) 

X 
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Figure 19: Axial-vector current correlators at am = 0.0016 (left) and 0.008 (right) for \Q\ 
= 2. The dashed lines are the qChPT prediction with parameters determined through \Q\ 
= and 1 sectors. 



{V{t)f = 2 X:(l + sYiP^x + Xo)P^(xo))l, (7.19) 

X 

at am = 0.0016, 0.0032, 0.0048, 0.0064, and 0.008. We take an average of the source point xo 
over {L/2Y x (T/2) lattice sites. In qChPT formulas ()7.3j) and ()7.4|1 . we have five parameters 
to be determined: F,r, S, Egg, m^ and a. Since these correlators are weakly depending on 
Ftt, we use the jackknife samples of obtained from the axial- vector current correlator, 
= 98.3(8.3) MeV. Unfortunately, there still remain too many parameters to fit with qChPT 
expressions. Therefore, we use the relation ()C13|) and an input TqX = ''^oiQ'^)/^ = 0.059(3) 
from a recent work [117], which gives mo = 940(80) (23) MeV, where the second error reflects 
the error of TqX- 

With these input values, we fit the correlators (j7.18|) and (j7.19|) in the range 7 < t < 13 
at different Q and m simultaneously. Fig. QUI shows the correlators with fitting curves. For 
\Q\ < 1 sectors, the data at all available quark masses am = 0.0016, 0.0032, 0.0064 and 
0.008 are fitted well, and we actually get x^/dof ~ 0.7. (Note that the correlations between 
different t's, m's and channels (PS and S) are not taken into account.) This fit yields S^/'^ = 
257 ± 14 ± 00 MeV, which is consistent with Ref. [96], S^^^ = 271 ± 12 ± 00 MeV, and a 
= —4.5 ± 1.2 ± 0.2, where the first error is the statistical error and the second one is from 
uncertainty of (Q^). 

Here are some remarks. The ratio 

^ = 1 + ^ (m^G(O) + «A(0)) , (7.20) 
should indicate the size of the NLO correction in the e expansion. Our data, 1.163(59), 
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implies that the e expansion is actually converging, a has large negative value, which is also 
reported in Refs. [118] and [35]. These results contradict with a previous precise calculation 
[119], which obtained a smaller value a = 0.03(3). If we instead assume a = 0, and fit 
the data with as a free parameter, we obtain F.^ = 136.9(5.3) MeV and S and Sefr are 
almost unchanged. (Detailed numbers are summarized in Table ^2) A possible cause is 
that \Q\ = 1 is not small enough to derive the partition function Eq. ()C.12|) (See Appendix 
Eq. ()C13|) may also have a systematic error due to finite V as well as finite a. Look at 
the data with \Q\ = 2, which are plotted in Fig. 1201 They do not agree with expectations 
from the qChPT shown by dashed curves in the plots, which is also seen in the case with 
the axial- vector correlators. A simultaneous fit with all the data including \Q\ = 0, 1 and 
2 gives a bad x^/dof (~ 12). Thus our numerical data might be posing a problem whether 
the partition function of qChPT Zq is reliable or not, when IQl/iQ"^) <^ 1 is not satisfied 
very much; (Q^) ~ 4.34(22) on our lattice. 
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Figure 20: Scalar (top) and pseudo-scalar (bottom) correlators at am = 0.0032 (left) and 
0.0064 (right). The dotted lines are the fit results with all available mass parameters am = 
0.0016, 0.0032, 0.0064 and 0.008 for topological charges and 1. 
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7.3.3 Chiral condensates 



Consider the low-mode contribution to the scalar and pseudo-scalar condensates with a fixed 
topological charge Q, 

= -(trZ}„i(x,x))« 



/ °" 1 \ \ 



:^ (1 — am/2) A j + m 
(7/'757/'(x))^ = -(tr75D„i(x,x)) 




Mow ]^ 



\ ( 1 — am/ 2) Aj + m 



i=l 




(7.22) 



If A^iow is large enough, the higher mode contribution, AD^^(a;,x), should be insensitive to 
the link variables U^iy) if \y — x\ ^ and thus to the global structure of the gauge field 
configuration, such as the topological charge. We, therefore, expect that the difference of 
the scalar condensates with different topology can be approximated with low-modes, 

- ((^^(a;))« - (V5^(a;))°) ~ - (#(^))?ow)- (7-23) 

In fact, Fig|2Il shows that the convergence up to A^iow = 200 + \Q\ is really good. For the 
pseudo-scalar condensate, the situation is much easier, since the condensate is determined 
by the zero modes only; 

- (^75^)q = J7. (7.24) 

Note that the contributions from non-zero eigen-modes cancel because of the orthogonality 
among different eigenvectors. As shown in Fig.|22l our data with iViow = 200-1- IQI low-modes 
perfectly agree with this theoretical expectation. 

The free parameter in the scalar condensate is Seg as seen in (17.11) and ()(I16jl . We 
compare our numerical result of Eq. ()7.23p with the qChPT result T,q{^') — Sq=o(/^')- We 
use the low-mode averaging over x (T/2) lattice sites. Fig. EHl shows the results as 

a function of quark mass for \Q\ = 1, 2 and 3. The lattice data agree remarkably well with 
the qChPT expectation with Scg = 271(12) MeV as determined from the (pseudo-) scalar 
connected correlators as presented in the previous section. In fact, if we fit these data with 
Seff as a free parameter, we obtain 256(14) MeV which is consistent with the result above. 
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Figure 22: Pseudo-scalar condensates approximated with A^iow = 200 + \Q\ low-modes. The 
lines represent the expectation ()7.24p . 
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0.0032, 0.0048, 0.0064, 0.008. The lines are qChPT predictions with S^^^ = 271(12) MeV. 



7.3.4 Disconnected PS correlators 



Let us rewrite the disconnected pseudo-scalar correlator in a fixed topological sector, 

/ \ Q 

(^(i)^(0))L = ( E {i5D;^\x, t; x, t)) tr (ts^^^'IO, 0; 0, 0)^ 



200+IQI 



1 



am/2)Xi + m 



Vi{x)vl{x) 



200+IQI 



xtr 



am/2)Xi + m 



+ 2 (tr (75AZ^„^ 
X ( tr 



200+IQI 



i=l 



+ 



tr 



am/2)Xi + m 

2 



t;,(0)^;|(0) 



(7.25) 



where x = {x,t). Here we assume that higher mode's contribution does not have correlation 
with any local operator 0{y) separated enough from x, i.e. 



{AD-\x,x)Oiy))' 



{AD;^\x,x)yx{0{y)Y 



(7.26) 



where the expectation value (■ ■ ■)' represents insensitivity to the topological charge. We also 
use the translational invariance {0{x)) = (0(0)). Since we know (tr (75AD~-'^(a;, x)))' = 0, 
the above equation should be well approximated with low-modes. In fact. Fig. |^ shows 
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Figure 24: Convergence of the pseudo-scalar disconnected correlators for one sample con- 
figuration at am = 0.008 for Q = (left) and Q = 1 (right). 



a very good convergence with the lowest 200+ IQI eigen-modes. Similar results were also 
obtained previously in the study of the rj' propagator with the Wilson fermion [111] and 
with the overlap fermion [120]. 

In the qChPT the correlator is written as 



(P'^(t))« = Jd^x{l + sy{2P''{x)P^{0)-P''{x)P%0)y 



C'P-^ f#KG(x) + aA{x)) - e+A(x) 



(7.27) 



where 



-Sv^^"*- -4 (1 + ?J) ^ -2 ((^) - ^) . (7,28) 

In Fig. I2S1 we present the data in \Q\ = 0-3 topological sectors at two representative 
quark masses am = 0.0032 and 0.0064. The qChPT predictions are plotted with the pa- 
rameters determined from the axial-vector and (pseudo-) scalar connected correlators: S^^^ 
= 257 MeV, = 98.3 MeV, mo = 940 MeV, and a = -4.5. One finds that the agreement 
is marginal, though the correlator's magnitude and shape are qualitatively well described. 
Instead, if we fit the disconnected correlator with S and a as free parameters while fixing 

and mo to the same value, we obtain T,^^^ = 227(32) MeV and a = —3.5(1.2), which are 
statistically consistent with those input numbers. We conclude that not only the connected 
correlators, but also disconnected correlators can be consistently expressed by the qChPT 
in the e-regime when \Q\ is small. Details of the fit results are listed in Table ITTl 
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Figure 25: Disconnected pseudo-scalar correlators in the < \Q\ < 3 sectors at am 
0.0032 and 0.0064. The curves represent the results of qChPT with E^/^ ^ 257 MeV, 
98.3 MeV, tuq = 940 MeV and a = -4.5. 



49 



Table 11: Summary of the fitting results in the e-regime. The first column denotes the topological sectors used in the fit. 
The values in [• • •] are input parameters. The first error is statistical. The second and third errors refiect the uncertainty 
in the input parameters, (Q^) and F^, respectively. 
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< IQI < 1 


125.7(5.6)(0.9) 


223(29)(00) 




[0] 


[734(33) (14)] 


0.7 


< IQI <2 


[98.3(8.3)] 


229(33)(00)(03) 




-3.6(0.2)(0.3)(1.0) 


[940(80) (23)] 


1.0 


0< IQI <2 


135.0(4.9)(1.4) 


237(32) (00) 




[0] 


[684(25) (13)] 


1.7 


< Q < 3 


[98.3(8.3)] 


229(33)(01)(05) 




-3.6(0.1)(0.2)(0.8) 


[940(80) (23)] 


1.0 


0< IQI <3 


139.3(4.1)(1.4) 


244(32) (00) 




[0] 


[663(19) (12)] 


1.9 



scalar condensate {{ipip)^ — (ipip)'^) 



1 < IQI < 3 256(14) L2 



8 Conclusions and discussions 



The exact chiral symmetry is established on the lattice, with the overlap Dirac operator which 
satisfies the Ginsparg- Wilson relation. However, the overlap Dirac operator becomes ill- 
defined at certain points where the zero-mode of Hw appears. Also, practically, these points 
make the numerical study very difficult. One has to carefully evaluate the discontinuity 
of the overlap fermion determinant (reffection/refraction) and the polynomial or rational 
expression of the overlap Dirac operator itself goes worse near this discontinuity. In fact, 
these dangerous points are topology boundaries; it is known that crossing Hy/ = changes 
the index of the overlap Dirac operator, or topological charge. 

In the continuum limit, = is automatically excluded but it would be better if we 
can construct the lattice gauge theory which does not allow Hw = from the beginning 
with a finite lattice spacing. There have been proposed two promising strategies. One is the 
gauge action which satisfies Liischer's "admissibility" condition and the other is an additional 
fermion determinant, det H^r. Both of them are designed for the use of the hybrid Monte 
Carlo algorithm, in which global and small (smooth) updates are performed. Therefore, 
we regard them as "topology conserving actions". In this thesis, we have investigated the 
possibility of lattice QCD in a fixed topological sector. We studied the properties of the 
topology conserving actions with no light quarks (namely, pure SU{3) theory.). Although 
the admissibility condition with small e (< 1/20) can strictly prohibit the topology changes, 
our interest is the case with e ~ 0(1) for practical purposes. In the (quenched) Hybrid Monte 
Carlo updates, we found that the topology change is strongly suppressed for 1/e = 2/3 and 
1, compared to the standard Wilson plaquette gauge action. The topological charge seems 
more stable for finer lattices, and it is possible to preserve the topological charge for O(IOO)- 
0(1,000) HMC trajectories at a ~ 0.08 fm and L ~ 1.3 fm, which might be applicable to 
the e- regime. The action Eq. ()2.15|) has, thus, been proved to be useful to accumulate gauge 
configurations in a fixed topological sector. While the gauge action Eq. (j2.15|) with 1/e = 1 
and 2/3 allows the topology changes, at times, the topology conservation with the fermion 
determinant detH^ seems perfect. The numerical cost for this determinant is, of course, 
much more expensive than the quenched case, but it would be negligible compared to the 
cost of the dynamical overlap fermion. The full QCD overlap fermion simulation with the 
additional determinant det Hy^ may be very efficient, if we can omit the refiection/refraction 
procedures. 

To test their practical feasibility, we measured the heavy quark potential with these ac- 
tions. The lattice spacing is determined from the Sommer scale Tq. With these measurements 
we also investigated the scaling violation at short and intermediate distances. The probe in 
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the short range is the violation of the rotational symmetry, and a ratio rc/ro of two different 
scale can be used for the intermediate distance. For both of these we found that the size 
of the scaling violation is comparable to the case with the Wilson plaquette gauge action, 
which is consistent with the expectation that the term with 1/e introduces a difference at 
most O(a^) or the additional large negative mass Wilson fermions are decoupled. These 
actions show no disadvantage as far as Wilson loops are concerned. We also found that 
the perturbative expansion of the coupling (after mean-field improvement) shows very good 
convergence even if 1/e term is introduced. The coupling constant in a certain scheme at a 
given scale is consistent among different values of 1/e. 

As another advantage of the (approximate) topology conservation, the low-lying eigen- 
values of the Wilson- Dirac operator in the negative mass regime is suppressed. This reduces 
the cost of the numerical implementation of the overlap-Dirac operator. We observed that 
the topology conserving actions have a gain about a factor 2-3 at the same lattice spacing 
compared to the standard Wilson gauge action. Comparison with the other improved gauge 
actions, such as the Liischer-Weisz, Iwasaki and DBW2 would be interesting. 

In a fixed topological sector, one of the very interesting applications would be the QCD 
in the e- regime. In the chiral hmit, the (q)ChPT analysis shows that the meson correlators 
are largely affected by the fermion zero-modes, and thus by the topology of the background 
gauge field. In order to study lattice QCD in such a regime, the chiral symmetric Dirac op- 
erator is essential, otherwise the large lattice artifacts would contaminate the fundamental 
points of the analysis, such as the definition of the topological charge, or what is the zero- 
mode, etc. To demonstrate how much effective the overlap Dirac operator is, we studied the 
quenched QCD with very small quark masses in the range 2.6-13MeV. In this chiral regime, 
we observed that the chiral behavior is really nice, as seen in, for example, accurate deter- 
mination of Za.1 or the pseudo-scalar condensates We also found the importance of 
the low lying modes and used them to extract the low-energy constants. (Note that the low 
lying modes are not affected by the a — > limit, since it would add higher modes only, which 
are irrelevant to the low energy dynamics.) Prom triplet meson correlators with Q = and 
1, we extracted = 98.3(8.3) MeV, = 257(14)(00) MeV i^'^ = 271(12)(00) MeV), 
mo = 940(80) (23) MeV, and a. — —4.5(1.2) (0.2). In these numerical results the second error 
reflects the error of (Q^)- We also obtained consistent results from disconnected pseudo- 
scalar correlator and the chiral condensates. Although we observed minor inconsistencies 
and problems due to finite Q correction, which is a special restriction of the quenched ChPT, 
these remarkable successes would encourage us to go further to A^/ 7^ full QCD in the e- 
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regime. 



The exact chiral symmetry with the overlap Dirac operator is surely a most remarkable 
progress in the lattice gauge theory and it would be more and more important in both of the 
numerical works and theoretical works, in the future. In order to perform the path integrals, 
one should, however, note the fact that there is several points where Hw has zero modes, 
which makes the overlap Dirac operator ill-defined, its locality is doubtful, and the numerical 
cost is suddenly enhanced. Our observation shows that the lattice QCD with fixed topology 
would be one of interesting and promising solutions to this. A number of applications such 
as 6'-vacuum, finite temperatures, and so on, might be possible as well. 

For the future works, we would like to give a few remarks. 

We have observed a large /3-shift when the negative mass Wilson fermion Eq. (j2.18|) is 
added, which may also cause a unwanted large scaling violation. We thus propose another 
"topology stabilizer" which would have effects on Hw ~ modes only; 

^2 2 ^ (8-1) 

where the denominator is corresponding to the twisted mass ghost [121-123] 

with a large negative mass — (1 + s)/a and a small twisted mass m^. The effect from both 

determinant should be canceled unless Hw has small eigenvalues. Note that det H^/^H"^ + 

a?mf) —>■ 1 with rrit fixed. 

As a final remark, let us consider the contribution from different topological sectors. 

In each topological sector, one measures the expectation value of operator O with a fixed 

topological charge, 

where dU^ denotes the integral over the gauge fields with topological charge Q and det{D + 
m)'^f is the determinant of Nf flavor overlap fermion with quark mass m. In order to 
calculate the total expectation value (in 6^- vacuum); 

, ^ Eq e'^Q / dU^O detjP + mffe-'o Eq e'^Q(0)QZQ/Zo 
EqC^'"^ JdUQdetiD + m)^fe-So Ege^^^^Q/^o ' 

the ratio Zq/Zq has to be evaluated. In fact, there are several proposals to calculate this 
ratio of the partition function, including the tempering method [124-128]. Here, we would 
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like to propose an easier method. With an assumption that the ratio should be Gaussian; 

Zq/Zo oc exp (- 2^^^g=o ) ' (8-4) 
in a large volume V, one only has to calculate the topological susceptibihty [129, 130], 

to evaluate the ratio Zq/Zq with any value of Q. In principle, the topological susceptibility 
can be evaluated with the local topological charge density operator q{x) — —Trj^aD{x, x)/2 
in a semi-local region (assuming the cluster decomposition principle); 

X = / d'x{q{x)qm'=' = j^^ d'x{q{x)q{0)y=', (8.6) 

where Jy, d^x is performed over the small volume V around the origin, where V ^ I/Kqcd- 
Thus, if the topological susceptibihty in a fixed topological sector, 

[ ^d'x{q{x)q{Q)f, (8.7) 

shows sufficiently small Q dependence and the finite volume effects due to V and V are 
both negligible, then should be a good approximation of % ^. What we would like 
to emphasize here is that summing up all the topology with the weight Zq/Zq (in a 9- 
vacuum) is not so difficult and maybe not so important unless 9 term is considered^, since 
Zq/Zq —>-v^oo 1, although the numerical studies has to be done in the future works to verify 
this quite optimistic argument. It seems an appropriate way to preserve the topology along 
the simulations, since any lattice gauge action would eventually create large barriers between 
the topological sectors in the continuum limit. 
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A Notations 

Lorentz indices /i, v run from to 3. The lattice spacing and size are denoted by a and L 
respectively. The gauge field U^j,{x) e SU{3) is located on the hnk from x to x + jla, where 
fi is the unit vector in direction The plaquette is denoted by Pnvix) — Un{x)Ui,{x + 
ap,)Uj^{x + av)Ul{x). The fermion field '^{x) is located on the site x. The forward and 
backward covariant difference of the fermion field ^^{x) are defined by 

V ^(x) = U^{x)i}} {x + jia) -ii){x) 



a 



, , i^ix) — UUx)ib(x — ua) 
Vlip{x) = ^ ^^Y^ ■ (A-l) 

Euclidean Dirac matrix 7^ satisfies 

li = {7m' Iv} = 25^,^, 75 = 7o7i7273- (A.2) 
The expectation values of an operator O with a fixed topology is denoted by (O)*^. 



B The hybrid Monte Carlo algorithm 

We explain how to generate the configurations in our numerical simulations. The hybrid 
Monte Carlo algorithm (HMC) is one of the most efficient [25] algorithms and widely used 
in the lattice QCD studies. 

Consider updating fields, g-j's, which have an action S{q). The HMC algorithm consists 
of two parts; 
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1. generate a candidate of updating field qf. 

2. judge whether qf is accepted or rejected. 

The former is called molecular dynamics steps, since it is very similar to the trajectory of a 
classical particle, randomly walking in configuration space, of which path is determined by 
the equation of motion; 

where the initial momentum 7rj(T — 0) is randomly given with probability density 



^({^.(0)}) = (n^exp(-X:^.'(0)/2)) . (B.2) 



The final point, qi{T — Tf), obtained in this way (after a fixed 'time' t/), is chosen to be 
the candidate qf. This step is known to satisfy "the detailed balance" which is a sufficient 
condition for an algorithm to sample the distribution with the correct Boltzmann weight 
ex.p{—S{q)). However, in the numerical studies, the evolution is done by the leap-frog updates 
with a finite step-size Ar, 

qi (r + At) - (t) = At * tt^ (r + Ar/2) , 
7ri(T + At/2) -7ri(T- At/2) = -AT*^^(r). (B.3) 

Therefore, one needs the latter procedure, or the Metropolis test, to judge that the new 
configuration is accepted or not; q!i is accepted with probability 

-7Tf/2-S{q'=) ^ 
^' g-7r|/2-5(g) 



P^min }, (B.4) 



which realizes the exact detailed-balance. The above steps {N^ds molecular dynamics steps 
with the step-size At and the Metropohs test) are done in every 'trajectory'. In this way, 
the configurations are generated performing the hundreds or thousands of trajectories, A^trj- 

Let us consider the path integral 

Z ^ JdU^ expi-Sc) det{D + 

= J dU^eM-Sa) J d(f)d(l)*exp j d'^x(f){x)*{{D + m)Y'^(f){x)^ , (B.5) 

where the pseudo-fermion field is introduced to evaluate the fermion determinant det{D + 
m)^, which can not be directly calculated when the lattice size is large. 
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Since the pseudo-fermion part of the action is a simple bihnear, can be updated with a 
conventional heat-bath algorithm. The HMC algorithm is applied to updating the link fields 
f/^'s with fixed cj) = (ph which is generated by the heat-bath method; 



Thus, our lattice QCD simulation is performed as follows, 

1. Choose a starting link variable configuration. 

2. Choose 7] = [D + m)^^(f)h to be a field of Gaussian noise (heat-bath). 

3. Choose the momentum of link variables, 7r^(x) from a Gaussian noise. 

4. Nmds molecular dynamics with the step-size Ar are done. 

5. Accept or reject the new configuration (Metropolis test). 

6. Store the (new or old) configuration. 

7. Return to step 2. 
Note that both of 



involves the calculation of (D + m)~^, which requires a very expensive computational cost 
in the chiral limit, m — > 0, since a lot of low eigenvalues appear. Also, it is notable that 
when {D + m) is not smooth with respect to the gauge link field ?7^, sudden jumps of 
momentum evolution in Eg. ()B.3jl should be monitored, otherwise over(under) estimation of 
the momentum would break the algorithm. Finally, we should refer to the fact that each step 
of HMC is small variation of the fields, and therefore, tends to keep the initial topological 
properties of the link variables. 

C Quenched chiral perturbation theory in the e-regime 

This appendix is devoted to a brief review of quenched chiral perturbation theory (qChPT) 
in the e-regime [107] and we summarize the relevant formulas for our analysis of meson 
correlation functions. 




(B.6) 



S'^'ff, and 




(B.7) 
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The partition function of qChPT with A^^ valence quarks is defined 

Z{e,M) = J dU exp J d'^xCliixij, (C.l) 
where the Lagrangian is given by 

= ^Strid^Uix)''d^Uix))-'^StT{UeUix) + Uixy'U,') 

9 

+^<f (^)' + ^d,^x)d,<!>{x), (C.2) 

at the leading order of pion mass, m^, and pion momentum, p^, expansion, where A''^ de- 
notes the number of color, F^^ is the pion decay constant and S is the chiral condensate. 
U{x) is integrated over a sub-manifold of the super-group Gl{Nv\N^), called the maximally 
symmetric Riemannian sub-manifold which is characterized by a matrix of form 

U=(^^ d)' ^ef/(iV,), DeGl{N,)/U{N,), (C.3) 

and Grassmannian A^^ x matrices B and C. Str denotes the super-trace. The quark mass 
is involved in the mass matrix M = {ml^ + ml^) where and ly are the identity matrix 
in the fermion-fermion and boson-boson blocks respectively. The 6 term which violates CP 
symmetry, enters through Ug = exp{i9/Ny)Iiy^ + 1^^. Note that in the quenched approx- 
imation the singlet field = ^Str[— z In ?7(a;)] is not decoupled, and its mass ml and 
coupling a are introduced [132]. 

The e-regime [43-48] is a special regime where the quark mass is so small that the 
pion Compton wavelength ~ I/wItt is larger than the linear extent of the space-time L. In 
this regime, the zero mode of pion gives important contribution and one must explicitly 
integrate out the constant mode of U{x). The other mode's contributions are perturbatively 
evaluated with the expansion parameter, ~ m-^/AnFT, ~ l/{LF^y. Namely, separating 
the zero mode and the other modes; 

f/(x) = f/oexp^^^SM, (C.4) 

the Lagrangian Eq. ()C.2|) can be rewritten. Together with e-expansion, one obtains the 
partition function with a fixed topological charge Q by Fourier transforming I^C.l^ as follows, 

Zq{m) = ^ r dee'''^z{e,M) 

= ^ f/O J dUl,daSdetU[,fexp - J d'x (c'^{x) + 
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1 

2^ 



X exp 



d9 J dU^d^iSdetUof exp 
'rnLV 



Vml ( FJ' 



V2iQ 



V2, 



-Str([/^ + f/, 



a 



1 



f 1 



niTy 



Str(t/^ + t/o 



a 



2N, 



where we use 



U 

$0 



Uoe 



-^StT{-ilnUo), 
UeUo, 

(Sdet?7^)Qexp 

^StrHln[/^) = <|.o+ ^. 



,4 V2iQ 



$0 



In the last hne of ()(I5j) . we perform 6 integral as a Gaussian; 



1 

— / dOexp 



exp 



27r(Q2) 



exp 



1 /■'^ 

— / (i^^' exp 

27r i-TT 



(C.5) 

(C.6) 
(C.7) 
(C.8) 
(C.9) 

(C.IO) 



(C.ll) 



where (Q^) = VmlF^/2N^ and 6' = 6 - V2%/F^. We need a condition, < 1, 

in order to justify this Gaussian integral, otherwise the integral Eq. ()G.ll|) should depend 
on $Q, which means that $o and 9 can not be treated independently and the last line of 
the partition function Eq. ()G.5|) is invalid. Therefore, all the results shown below are reliable 
only for small \Q\. 

With redefinition Uq = Uq, the partition function in a fixed topological sector is obtained 

Zq{M) = ^ r^d9e'''^Zie,M) 

zn J-TT 
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27r(g2) J '12' 



(C.12) 



where dUo denotes the Haar measure of the maximally Riemannian sub-manifold of 
Gl{Ny\Ny). Note that the topological charge distributes as a Gaussian with variance 

which is in a good contrast with the full theory, for which (Q^) = mT,V/Nf is expected with 
Nf flavors. 

In the following we consider = 1 and 2 cases only, as wc arc interested in the system 
with two light quarks. Any correlators or condensates arc obtained by the perturbation of 
^ fields and the exact integration over zero mode, Uq, which can be written in terms of the 
Bessel functions. 

First let us calculate the condensates. At the tree-level the scalar condensate is given by 

- {{^^^^;)Q = Egipi) = E/.(7|Q|(/.)X|Q|(/.) + 7|q|+i(a^)X|q|_i(/x)) + (c.14) 

A'' 

with ji = mTlV. IiQii/j) and K\Q\{fj,) denote the modified Bessel functions. One- loop contri- 
bution does not change its functional form [133] 

E^-^°"P(/.) = Eeff/.'(/|Q|(/.')i^|Q|(/^') + + Seff^ = (C.15) 

A* 

but the parameters //. and S are shifted to ii' and Eefj: 

=mEefrl^, Eefr = E(H I. (C.16) 

Here, parameters (5(0) and A(0) are ultraviolet divergent tadpole integrals, 

which need renormalization. 

Next, we consider the fiavor-singlet meson operators 

S%x) = i,{x)Ii,y^{x), P\x) = i,{x)i-i^lN^^P{x). (C.18) 
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To add these operators to the QCD Lagrangian as source terms; 

£ ^ £ + s{x)S'^{x) + p{x)P%x), 
corresponds to the substitution 

M ^ M + s{x)In, + ip{x)lN^ 



(C.19) 



(C.20) 



in qChPT. The two-point correlation functions are obtained by differentiating the generating 
functional by s{x) and p{x). To O(e^) the results are 



^2 r 



{s\x)s'm 



Q 



where 



2F2 



Nr. 



{mlG{x) + a/\{x)) - A{x) 
{mlG{x) + aA{x)) - A{x) 



a 1 + a_ — 4 



a+ + a_ + 4' 



+ 1 + ^ 



and the constant terms arc given by 



c% 



y2 

^cff ^1-loop _ v2 



a_ =4 



// E ^ ^ -2 



(C.21) 
(C.22) 

(C.23) 



+ 1 + 



^eflf 1-loop v^2 



12 



/2 



1 Sq(^o 



/2 



(C.24) 
(C.25) 



Note that the prime denotes the derivative with respect to /i, 

= IMKM - I\Q\MK\Q\-i{n) - ^. (C.26) 

For flavor non-singlct case, we need a A^^ = 2 super-group integral, which is also described 
with the Bessel functions. The non-singlet operators are given by 



P%x) = '4^ix){Ty2)i^5lNyipix), 
with the Pauli matrices r". To O(e^) the two-point functions are given by 



{S'\x)S'\0)) 



Q 



— {mlG{x) + aA{x)) - A{x)b. 



{P'^{x)P%0))q ^ C%- 



2F2 



^{mlG{x) aA{x)) - A(x)6-F 



(C.27) 
(C.28) 



(C.29) 
(C.30) 
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where 



2 1 + 



c_ = -2 



(C.31) 



and 



2 \ Seflf / 

For the flavor non-singlet axial- vector current 



eflf 



yl«(a;)=^(a;)(r72)^7^75^(a;), 



the correlator is obtained [108] 



(A«(x)>1«(0))q = __^_2mEQ(//)A(a;). 



(C.32) 



(C.33) 



(C.34) 



An important observation here is that the axial-current correlator does not depend on the 
parameters related to the quenched artifact, i.e. rn^ and a. Also note that the constant term 
is proportional to rather than S. Therefore, this channel is suitable for an extraction of 
Ftt, whereas for the pseudo-scalar and scalar correlators appears only in the coefficients 
of A{x) and G{x) terms. 

For zero spatial momentum projection, it is convenient to define 



r2 1 T i 30 



(C.35) 
(C.36) 



As a Remarkable feature in the e-rcgime, the correlators do not behave as the usual expo- 
nential fall-off exp(— Mt) with the mass gap M. 
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